Fermion self-energies and pole masses at two-loop order in a general renormalizable 

theory with massless gauge bosons 



Stephen P. Martin 
Physics Department, Northern Illinois University, DeKalb IL 60115 USA 
and Fermi National Accelerator Laboratory, PO Box 500, Batavia IL 60510 

I present the two-loop self-energy functions and pole masses for fermions in an arbitrary renormal- 
izable field theory, in the approximation that vector bosons are treated as massless. The calculations 
are done simultaneously in the mass-independent MS, DR, and DR renormalization schemes, with 
a general covariant gauge fixing, and treating Majorana and Dirac fermions in a unified way. As 
examples, I discuss the two-loop strong interaction corrections to the gluino, neutralino, chargino, 
and quark pole masses in minimal supersymmetry. All other two-loop contributions to the fermion 
pole masses in softly-broken supersymmetry can also be obtained as special cases of the results given 
here, neglecting only the electroweak symmetry breaking scale compared to larger mass scales in 
two-loop diagrams that involve W or Z bosons. 
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I. INTRODUCTION 

The CERN Large Hadron Collider and a future 
electron-positron linear collider should discover and, to- 



gether, thoroughly explore [1] the mechanism behind 
electroweak symmetry breaking. The small ratio of the 
scale of electroweak symmetry breaking to the Planck 
mass scale suggests that supersymmetric particles will 
also be found at these next-generation experiments. If 
so, then a primary goal of both experimental and theo- 
retical research will be to unravel the mechanism behind 
supersymmetry breaking. The most important clues will 
be the masses of the superpartners and the Higgs scalar 
bosons. Therefore it is important to be able to compute 
the physical masses accurately in terms of the underly- 
ing Lagrangian parameters, including at least the leading 
two-loop effects. 

In this paper, I will present results for the two-loop 
contributions to fermion self-energy functions and physi- 
cal pole masses in a general renormalizable field theory, in 
terms of the running renormalizcd couplings and masses. 
The approach used is intended to be as flexible as pos- 
sible, so that a common framework of calculation can 
be used to treat both Majorana and Dirac fermions, in- 
cluding chiral interactions, in both supersymmetric and 
non-supersymmetric theories. As a simplifying approxi- 
mation, vector bosons will be treated as massless in the 
two-loop parts in this paper. In the Standard Model and 
extensions of it that do not enlarge the gauge group, this 
amounts to neglecting the effects of electroweak symme- 
try breaking compared to the masses of heavier particles 
in two-loop diagrams that have W and/or Z boson prop- 
agators. (The effects of non-zero W and Z boson masses 
can be included as usual in the one-loop part.) This will 
likely be a good approximation for the pole masses of 
the top quark and most of the supersymmetric particles, 
because of the exclusions of light squarks, sleptons, and 
gluinos already achieved by the CERN LEP e + e~ collider 
[2] and the Fermilab Tevatron pp collider [3, 4] . 

The mass defined by the position of the complex pole 
in the propagator is a gauge-invariant and renormaliza- 
tion scale-invariant quantity [5]-[13]. The pole mass in 
principle does suffer from ambiguities [14] due to infrared 
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physics associated with the QCD confinement scale, but 
these are probably relatively too small to cause a practi- 
cal problem for strongly-interacting supcrpartners. The 
pole mass should be closely related in a calculable way to 
the kinematic observable mass reported by experiments 
[15]. In recent years, many important higher-order cal- 
culations of self-energy functions and pole masses in the 
Standard Model have been performed, including two-loop 
[16]-[19] and three-loop [20]-[21] contributions for quarks, 
and two-loop results for electroweak vector bosons [22]- 
[28]. In addition, there are important two-loop results 
for top and bottom quarks [29-31] and the gluino [32] 
in low-energy supersymmetry. The general treatment of 
the present paper will confirm and extend the results of 
those papers. 

The notation and strategy used here are very similar to 
those found in my previous papers on scalar self-energy 
functions and pole masses at two-loop order in a general 
theory [33, 34] . In the next section, I review the conven- 
tions used, discuss the formalism for self-energy functions 
and pole masses for fermions in a two-component nota- 
tion, and review the methods used for numerical evalua- 
tion of the required two- loop integrals. 



II. NOTATIONS AND SETUP 

A. Notations for fields, interactions, and indices 

In this paper, the spacetime metric tensor is 
diag(-l,+l,+l,+l). 



n 



(2.1) 



I use a two-component notation for fermions, as in 
ref. [35] and similar to that found in ref. [36] . Left-handed 
spinor fields tp a always carry undottcd spinor indices 
a, (3,. . . = 1,2, and right-handed spinor fields always 
carry daggers and dotted spinor indices a, $, . .. = 1,2, 
with 



v>! = wo*- 



(2.2) 



However, the spinor indices are most often suppressed, 
as described below. The spinor indices are raised and 
lowered with the two-index antisymmetric symbol with 
components e 12 = — e 21 = 621 = —£12 = 1, and the 
same set of sign conventions for the corresponding dotted 
spinor indices. Thus 



"0a = £ Q /3'0 /3 , tp C 



c af3 



(2.3) 
(2.4) 



Spinor bilinears can be combined to form vector quanti- 
ties using the matrices a ^ a p and ~off defined by 



(Tq = Co 
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1 
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When constructing Lorentz tensors from fermion fields, 
the heights of spinor indices must be consistent in the 
sense that lowered indices must only be contracted with 
raised indices. As a convention, indices contracted like 



and a a , can be suppressed. For example, 



(2.6) 
(2.7) 
(2.8) 
(2.9) 



The behavior of the spinor products under hermitian con- 
jugation (for quantum field operators) or complex conju- 
gation (for classical fields) is as follows: 



£x = 

















(^X^=X*^, 



(2.10) 
(2.11) 
(2.12) 



The following identities also hold: 

[a^a" + a v a»]J = -2^6%, (2.13) 
[cfV + 0V] d / 3 = -2^£|, (2.14) 

Tr[a^a v ] = Tv^a"] = -2rf v . (2.15) 

In terms of two-component fermion notation, a single 
Dirac fermion is given in the chiral representation by 



td ' 



(2.16) 



where £ is the two-component fermion that describes 
the left-handed part of ^ and x is t ne two-component 
fermion that describes the conjugate of the right-handed 
part of ^. The Dirac matrices are 



7^ 



P 
Ju 



(2.17) 



In this paper, I consider a general renormalizablc 
field theory, containing t a set of real scalars Ri, two- 
component Weyl fermions ipi, and vector bosons V£. 
Scalar field indices are i, j, k, . . ., fermion flavor indices 
are /, J,K, . . ., and a, b, c, . . . run over the adjoint rep- 
resentation of the gauge group, while (i, v, . . . are space- 
time vector indices. Repeated indices of all types are 
summed over unless otherwise noted. 

The masses and couplings are evaluated by taking the 
fields in the Lagrangian in a squared-mass eigenstate ba- 
sis, after the Higgs fields are assumed to have been ex- 
panded around their vacuum expectation values as deter- 
mined by the loop-corrected effective potential (so that 



t A complex scalar can be written as two real scalars, and a Dirac 
fermion as two Weyl fermions, so this entails no loss of generality. 
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tadpole graphs do not contribute). The kinetic part of 
the renormalized tree-level Lagrangian is then written as: 



-ld u Rid»Ri - \m 2 R 2 



-i^o^d^i - -(m'-fy/^J + c.c.) 

-^(W - - \m 2 a V^. (2.18) 

The non-gauge interactions of the scalar and fermion 
fields are given by the renormalized Lagrangian: 



-—\ l i k RiRjRk — —\ l -' km RiRjRkRm 



-\{y JKi M K Ri + c.c). 



(2.19) 



where A yfe and A yfem are real couplings and the Yukawa 
couplings y JKl are symmetric complex matrices on the 
indices J, K , for each i. Raising or lowering of fermion 
indices implies complex conjugation of the Lagrangian 
parameters, so 



m u = (m 1J ) 



y JKl ee (y JK r 



(2.20) 



Actually, without loss of generality, m IJ can be taken to 
have only real and non-negative entries, but the index 
height convention is maintained for clarity. The heights 
of real scalar and vector indices have no significance, and 
in any given equation are chosen for convenience. 

The scalar squared masses mf and the fermion squared 
masses mjKm KJ — m 2 5j are taken to have been diag- 
onalized (by an appropriate rotation of the fields if nec- 
essary). However, the fermion mass matrix m IJ is not 
necessarily diagonal; instead it must have non-zero en- 
tries only when I and J label two-component fermions 
with the same squared mass and in conjugate represen- 
tations of the unbroken gauge group. In particular, when 
dealing with Dirac fermions, it is most useful to work in 
a basis in which the corresponding matrix m IJ contains 

m 



2x2 blocks of the form 







on the diagonal. 



Next consider the gauge interactions of the theory. Let 
T a be the Hermitian generator matrices of the gauge 
group for a (possibly reducible) representation R. They 
are labeled by an adjoint representation index a corre- 
sponding to the vector bosons of the theory, . They 
satisfy [T a , T b ] = if abc T c , where f abc are the totally anti- 
symmetric structure constants of the gauge group. Then 
results below are written in terms of the invariants: 



(T a T a ) l j 
Tr[T a T b ] 

jacd jbed 



c(i)6{, 

I{R)S ab , 
C(G)S ab , 



(2.21) 
(2.22) 
(2.23) 



summed over the representation R, and the Casimir in- 
variant of the adjoint representation of the group, re- 
spectively. When the gauge group contains several sim- 
ple or U(l) factors with distinct gauge couplings g a , 
the corresponding invariants are written C a (i), I a (R), 
and C a (G). The normalization is such that for SU(N), 
C(G) = N and each fundamental representation has 
C(i) = (N 2 - 1)/2N and contributes 1/2 to I(R) for 
each Weyl fermion or complex scalar. For a U(l) gauge 
group, C(G) = and a representation with charge q has 
C(i) = q 2 and contributes q 2 to I(R). The two-loop re- 
sults given below will be presented in terms of these group 
theory invariants for the representations carried by the 
scalar and fermion degrees of freedom. 

The preceding paragraph applies to the two-loop parts, 
in which the gauge group is treated as unbroken and 
m 2 = 0. In the one-loop parts of the self-energy func- 
tions and the fermion pole masses, the effects of non-zero 
vector boson masses will be retained. This means that 
the gauge group cannot be treated as unbroken, and the 
interactions of the vector bosons with the fermions have 
a more general form. They can be written as: 



-gauge 



-gfV^a^j 



(2.24) 



where g^ J are couplings obtained by going to the tree- 
level mass eigenstate basis for the fermions and vector 
bosons. In the special case of an unbroken gauge sym- 
metry, one has gf T = .g a [r a ] 7 J . 

The computations in this paper are performed with 
a vector boson propagator obtained by covariant gauge 
fixing in the usual way: 

-i5 ab ( V ^ + k»k»C ml ) [l/(fc 2 + m 2 a )] , (2.25) 

where for later convenience I use the notation 

C x f(x) = [f{x)-f(tx)]/x, (2.26) 
with the appropriate limit for massless vectors: 

Iim[£ x /(x)] = (l-0/'(0). (2.27) 

Here £ = 0, 1, and 3 correspond to the Landau, Feynman, 
and Fried- Yennie gauge-fixing choices, respectively. The 
self-energy functions depend on £, but the pole masses 
do not. For the two-loop computations below, the vector 
bosons are treated as massless, so the propagators are 



(2.28) 



Infrared divergences are dealt with by first computing 
with a finite vector boson mass, and later taking the 
massless vector limit. All contributions involving gauge 
boson loops implicitly include the corresponding contri- 
butions of ghost loops. 



B. Regularization and renormalization 



which define the quadratic Casimir invariant for the rep- 
resentation carrying the index i, the total Dynkin index 



For each Feynman diagram, the integrations over in- 
ternal momenta are regulated by continuing to d = 4 — 2e 



4 



dimensions, according to 



devices. First, 



/ 



d 4 k -» (27r M ) 2e / d d k. 



(2.29) 



In the dimensional regularization scheme, the vector 
bosons also have d components, while in the dimensional 
reduction scheme they have d ordinary components and 
2e additional components known as epsilon scalars. This 
means that the 4-dimensional metric in the vector prop- 
agator of eq. (2.25) is replaced by 



rT/(k 2 +ml) 



g^/(k 2 + m 2 a ) 

+g^/(k 2 + m 2 a + m 2 ), (2.30) 



where g^ is projected onto a formal d-dimensional sub- 
space, and <7 M1/ onto the complementary 2e-dimensional 
subspace, and m 2 is the epsilon scalar squared mass pa- 
rameter. (In general, there should be a different m 2 for 
each a, but it should cause no confusion to omit the ad- 
ditional subscript in this paper.) Counterterms for the 
one-loop sub-divergences and the remaining two-loop di- 
vergences are added, according to the rules of minimal 
subtraction, to give finite results, which then depend on 
the renormalization scale Q given by 



Q 2 = 4ire 



v. 



(2.31) 



Logarithms of dimensionful quantities are always written 
in terms of 



InX = \n(X/Q 2 ). 



(2.32) 



The resulting renormalization schemes are known as MS 
[37] and DR [38] , respectively, for the cases in which g^ v 
is not and is included. 

The epsilon-scalar squared mass parameter m 2 appear- 
ing in the DR scheme is unphysical. One could set m 2 
equal to zero at any fixed renormalization scale, but 
then it will be non-zero at other renormalization scales, 
since it has a non-homogeneous beta function [39] . Fur- 
thermore, under renormalization group evolution it will 
feed into the ordinary scalar squared masses in the DR 
scheme. Fortunately, a redefinition (given in [40] at one- 
loop order, and at two- loop order in [41]) of the ordi- 
nary scalar squared masses completely removes the de- 
pendence on the unphysical epsilon scalar squared mass 
m 2 from the renormalization group equations and the 
equations relating tree-level parameters to physical ob- 
servables in softly-broken supersymmetric theories. The 
resulting DR scheme [40] is therefore an appropriate one 
for realistic models based on supcrsymmctry, such as the 
Minimal Supersymmetric Standard Model (MSSM). In 
this paper, calculations will be presented simultaneously 
in all three schemes, using the following two notational 



^MS — 



(2.33) 



1 for MS 
^0 for DR, DR'. 
Second, terms that involve the unphysical parameter m 2 
should be construed below to apply only to the DR 
scheme, not the DR or MS schemes. 



C. Self-energy functions and pole masses for 
fermions using two-component notation 

The full, loop-corrected Feynman propagators with 
four-momentum p^ are denoted as shown in Fig. 1, which 
defines C/, D IJ , and Djj as functions of the masses and 
couplings of the theory and of the external momentum 
invariant 



s = —p 2 . 



They are given, starting from tree level, as 



Dij 

a 



= m 



i, j 



/(p 2 + m 2 ) + 
mu/(p 2 +mj) + 
5i/{p 2 + m 2 ) + .. 



(2.34) 



(2.35) 
(2.36) 
(2.37) 



with no sum on I in each case. In general, D IJ is a 
complex symmetric matrix, and Djj is obtained from 
it by taking the complex conjugate of all Lagrangian 
parameters appearing in its calculation, but not taking 
the complex conjugates of loop integral functions, whose 
imaginary (absorptive) parts correspond to fcrmion de- 
cay widths to multi-particle intermediate states. 

The computation of the full propagators can be or- 
ganized, as usual in quantum field theory, in terms of 
one-particle irreducible self-energy functions. These are 
defined in Fig. 2. (The same remark applies for the re- 
lationship between the functions Q IJ , Sl/j as for D IJ , 
Dij.) Then one has the matrix diagrammatic identity 
shown in fig. 3. To write this in terms of the self- 
energy functions, denote N x N matrices (where N is the 
number of two-component left-handed fcrmion degrees of 
freedom, so that /, J = 1, 2, . . . N): 



c/ = (c T yj 

B IJ = D IJ , 
£/ = (S T ) J 7 ee £/, 

n IJ ee tt 7J , n 



D 



D 



i.j, 



i, j 



n 



m 



m 



u 



m/j EE 771/ j, 



(2.38) 
(2.39) 
(2.40) 
(2.41) 
(2.42) 



Then fig. 3 implies that the propagator functions obey 
the 4N x 4N matrix equation: 
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— ► — 

J 




— ► 

I 



a 
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I 




— < — 

J 



lp 



a Al3 Cj 
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— ► — 

I 




— < — 

J 



a 
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■* — 

I 




— ► 

J 



iSj Du 



FIG. 1: The full loop-corrected propagators for fermions in two-component notation are associated with functions C/(s), 
D IJ (s), and Djj(s), as shown. Here s = — p 2 is the external momentum invariant. The shaded boxes represent the sum of all 
connected Feynman diagrams, and the external legs are included. 



V P 
► ► 




ip-a^zj/s ip-a^j/s -isJn IJ -iS^pSiu 

FIG. 2: The self-energies for fermions in two-component notation are associated with functions Ej(s) (for chirality-preserving 
propagation), and ft IJ (s) and flij(s) (for chirality- violating propagation), defined as shown. Here s = — p 2 is the external 
momentum invariant. The shaded circles represent the sum of all one-particle irreducible, connected Feynman diagrams, and 
the external legs are not included. 




FIG. 3: The diagrammatic version of the fermion self-energy identity in eq. (2.43). 




'ip-o [1 - S T /s] i[m + 0] N 
i[m+n] ip-a[l-'S/s] 



(2.43) 



The pole mass can be found most easily by considering the rest frame of the fermion, in which the space components of 
the external momentum p M vanish. This reduces the spinor-index dependence to a triviality. It follows from eq. (2.43) 
that the (complex, if the fermion is unstable) poles of the full propagator are the solutions for s of the non-linear 
N x N matrix eigenvalue equation: 

Det[sl-(l-S/.s)- 1 (m + n)(l-S T /s)- 1 (m + rj)] =0. (2.44) 

This can be solved iteratively by first expanding each of the self-energy functions in a Taylor series in s about the 
tree-level squared masses mj. Write the one- and two-loop contributions to the self-energy functions as: 

« " i^< ,J + a^ !,J + - (2 - 45 > 
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1 



an + . 



Hi--' ' 7 + (16tt 2 ) 2 " 



(2.47) 



where the superscripts (1) and (2) refer to the one- and two- loop contributions respectively. Then define the quantities 
(with sums on /', J', K, and K', but not on / or J): 



(i)j 



(2)J 



EfVj/i + mi7/S^ )J, m JJ 7s + m ir to {2)JI ' + il?},m JJ ' 



(2.48) 



+ {sf )Jf aWm JJ ' + fVfffi^' + m„,sffiW JJf + fi&E™ V }/s, (2.49) 
which play a role analogous to the self-energy functions of scalar or vector bosons, with eq. (2.44) taking the form 

sS J - (m 2 S J + — !— TT (1)J 4- 1 n (2)J Y 
Sdl V"'*' + 167T2 11 ' +(16^)2^ ). 



Det 



0. 



(2.50) 



It follows that the pole squared masses for the fermions are given by (with no sum on the index I): 



Sp i ei / = Mf - if/M/ = mj + 



16tt 2 



(16^ 2 ) 2 



n? )/ + n «' + E n<^nf/( m ? - m » )] , (2.5i) 



where one must put s 



ie (note with an infinites- 



imal positive imaginary part; this is necessary to give 
the correct negative imaginary part to the pole mass) 
everywhere on the right-hand side. Terms that are of 
three-loop order have been dropped. 

In writing eq. (2.51), it is assumed that the fermions 
that mix with each other are not degenerate, so that the 
last term is part of a well-defined perturbative expansion. 
If (nearly) degenerate fermions do mix, then the appro- 
priate version of (nearly) degenerate perturbation theory 
should be used instead to solve eq. (2.50). One can also 
obtain a solution iteratively, by first taking s — mj as the 
argument of the self-energy functions in eq. (2.50), solv- 
ing for s to obtain the next value for the argument Re(s) 
of the self-energy functions, and repeating until sufficient 
numerical convergence is obtained. However, despite the 
formal gauge invariance of the pole mass, this iterated 
procedure does not give a gauge-invariant result at two- 
loop order when massless gauge bosons are present, be- 
cause of the branch cut in the one-loop self-energy that 
is present except in the Fried- Yennie gauge £ = 3. This 
is because the pole mass result obtained by the iterative 
procedure is not formally analytic in the gauge coupling 
for £ 7^ 3, as explained in more detail in the analogous 
case for scalars in ref. [34]. 

For taking the limit s — > mj in eq. (2.51), it is conve- 
nient to define (again with no sum on I): 



ny 



(2) 



lim 



nf )J + 



mi 



ds 



(2.52) 



since this combination is independent of the gauge-fixing 
parameter £, and free of logarithmic divergences of the 



form ln(l — s/mj) that do appear in the individual terms 
when there are massless gauge bosons. The results for 
one-loop self-energy functions and pole squared mass 
contributions S^ J , fl^ IJ , flf), and Il[ 1 - )J , will be re- 
viewed in section III. The two-loop contributions to 
E (2)J q(2)ij (2) ; and fj(2) are presented in sec tion IV. 



D. The Feynman diagrams 

The one-loop and two-loop Feynman diagrams needed 
for the results just mentioned are shown in fig. 4. They 
are labeled according to a system described in ref. [33]. 



E. Two-loop basis integrals 

The results below will be written in terms of two- 
loop integral basis functions, following the notation given 
in [42, 43]. The one-loop and two- loop integral func- 
tions are reduced using Tarasov's algorithm [44, 45] to a 
set of basis integrals A(x), B(x,y), I(x,y,z), S(x,y,z), 
T(x,y,z), U(x,y,z,u), and M(x,y,z,u,v), correspond- 
ing to the Feynman diagram topologies shown in fig. 5. 
Here x,y,z,u,v are squared mass arguments. The addi- 
tional arguments s and Q 2 are not shown explicitly, be- 
cause they are the same for all functions in a given equa- 
tion. The functions A(x) and I(x, y, z) do not depend on 
the external momentum at all, with A(x) = x(lnx — 1) 
and I(x, y, z) = S(x, y, z)\ s =q. Each of the basis integral 
functions contains counterterms that render them ultra- 
violet finite. The precise definitions, and the calculation 
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Vfvvff Vfvvss Yfvvs 



FIG. 4: The one-loop and two-loop Feynman diagrams for fermion self-energies in the approximation of this paper. Dashed 
lines stand for scalars, solid lines for fermions, and wavy lines for massless vector bosons. Diagrams involving vector boson loops 
also include the corresponding ghost loop diagrams. The label for each diagram refers to a corresponding function obtained as 
the result of the two-loop integration. All counterterm diagrams for each diagram are included in these functions, rendering 
them ultraviolet finite. For each diagram, fermion mass insertions (indicated by adding a bar to the corresponding subscript F 
in the name) are to be made in all possible ways. Diagrams indicated by (*) vanish identically in the MS scheme with massless 
vector bosons, but not in the DR scheme with non-zero epsilon-scalar masses. 




A(x) B(x,y) I(x,y,z) S(x,y,z) T(x,y,z) U(x,y,z,u) M{x,y,z,u,v) 



FIG. 5: Feynman diagram topologies for the one- and two-loop self-energy basis integrals used in this paper. The letters 
x,y,z,u,v refer to the squared masses of the corresponding propagators. The dot on the T diagram means that T(x,y,z) = 
— dS(x, y, z)/dx. The precise definitions of these Euclideanized scalar integral functions, and methods for their evaluation, are 
described in [42, 43]. 
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of these functions and a publicly available computer code 
(TSIL) for that purpose, are described in [42, 43]. 

Several shorthand notations will be used. As explained 
in refs. [42, 43], it is convenient to define: 

T{0,y,z) = ]hn[T(x,y,z) + B(y,z)tox], (2.53) 
V(x,y,z,u) = -dU(x,y,z,u)/dy. (2.54) 

A prime on a squared-mass argument of a function is used 
to denote a derivative with respect to that argument, so: 

B{x,y') = dB(x iy )/dy, (2.55) 
I(x',y,z) = dl(x,y,z)/dx. (2.56) 

A prime on a function itself indicates a derivative with 
respect to the external momentum invariant s, so: 

B'(x,y) = dB(x,y)/ds (2.57) 
U'(x,y,z,u) = dU(x,y,z,u)/ds. (2.58) 

I 



But, note that below primes also appear on fcrmion in- 
dices, where they are used for a completely different pur- 
pose; fermions are labeled by indices / and I' if they 
combine to have a common squared mass rnj. 

Each of the functions in cqs. (2.53)-(2.58) can be 
reduced to combinations of other basis functions; see 
eqs. (3.1), (3.22), (4.14), (4.26), (5.3), and (6.18) of 
ref. [42] for formulas in the notation of the present pa- 
per. However, this explicit reduction is not done below 
in cases where it would needlessly complicate the expres- 
sions. 



III. FERMION SELF-ENERGY FUNCTIONS 
AND POLE MASSES AT ONE-LOOP ORDER 

In this section, I review the results at one-loop order. 
The chirality-preserving and chirality-violating parts of 
the fcrmion self-energy function are respectively: 



Z? )J - y IKi y JKi B FS (m 2 K ,m 2 ) + gi K g a K J B FV (m 2 K ,ml), (3-1) 

n WU = /«/^'« m ^B^(™ 2 K ,m2)-.g^m^'B 1?v (m^,m2), (3.2) 

where 

B FS (x,y) = [(y-x-s)B(x,y)-A(x)+A(y)]/2, (3.3) 

B Ts {x,y) = -B(x,y), (3.4) 
B FV (x,v) = (v - x - s)B{x,v) + A(v) - A{x) + sS m + £ v [{v(s + x) - (x - s) 2 }B(x,v) + (x - s)A(v)]/2, (3.5) 

B FV (x, v) = 3B(x,v)+ZB(x,Zv)-25 m . (3.6) 



These follow from direct evaluation of the first two Fcynman diagrams, with and without mass insertions, in fig. 4. [The 
result for Q^ 1 ] follows from cq. (3.2) by replacing the coupling parameters by their complex conjugates.] Here I have 
allowed for the possibility of general fermion-fermion- vector interactions and vector masses arising from spontaneous 
breaking of gauge symmetries. In the following, I will also make use of: 

B' FS (x,y) = [(y-x-s)B'(x,y)-B(x,y)]/2, (3.7) 
BL s ( x ,y) = -B'(x,y), (3.8) 

where the prime means a derivative with respect to s. 

In the special case of massless vectors corresponding to unbroken gauge symmetries, one makes the simplifications:^ 

9i K 9 a K J B F v(m 2 K ,m 2 a ) - g 2 a C a (I)SjB FV (mj,0), (3.9) 

-jtf™ M B fy (m| inl J) -> g 2 a C a {I)m u B Tv {m 2 ,Q), (3.10) 

where 

B FV (x,0) = £[s-{x + s)B(0,x)-A(x)}-s + s8 m , (3.11) 

B Tv (x,0) = (3 + OB(0,x)-26 m . (3.12) 



t The minus sign in eq. (3.10) occurs because the left-handed 
fermions with labels K and K' necessarily occur in conjugate 
representations of the unbroken gauge group. 
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It follows that the quantity defined in eq. (2.48) is 

n (i).7 = n (i,o).7 + n (i,iy (313) 

where the contribution from scalar exchange is: 

n^ 0)J = (yiKiy JKi m 2 J + m ir y I ' Ki yj, Ki m JJ ')B FS (m 2 K ,m 2 i )/s 

+ (m II ,y I ' Ki y JK ' l m KK , + y IKl yj> KHm KK ' 'm JJ ' ')B FS (m 2 K , mj), (3.14) 

and the contribution from vector exchange is: 

n (l,l)J = (gl K g« K > m Z+ miI ,gfg«/)B FV ( m 2 K ,m 2 a )/ S 

-(rrurgfgtfm™' + gi K g a /' m KK/ m JJ ')B FV (m 2 K ,m 2 a ). (3.15) 
In the special case of massless vector bosons, the latter expression reduces to: 

Ilj M)J = 2g 2 a C a (I)Sfmj[B FV (mj,0)/s + B FV (mj,0)], (3.16) 

with the well-known limit: 

Urn n< 1,1)J = 2g 2 a C a (I)Sjmj[5 -S m - 3Wf]. (3.17) 

s— >mj 

The above expressions can be inserted in the formula eq. (2.51) to obtain the one-loop contribution (and part of the 
two-loop contribution) to the pole squared mass. 



IV. FERMION SELF-ENERGY FUNCTIONS AND POLE MASSES AT TWO-LOOP ORDER 

In this section I present the results for the two-loop contributions to the self-energy functions and pole squared 
masses of fermions as defined in fig. 2. The results are divided into parts due to diagrams with no vector propagators, 
one vector propagator, and two vector propagators, with superscripts (2, 0), (2, 1) and (2, 2) respectively: 

s (2)7 = s (2,0)7 + s (2 > l)7 + s (2 > 2)7 5 (41) 
n (2)IJ = Q(2,0)JJ + Q(2,l)/J + fi (2,2)IJ (42) 
g(2) = fi (2,0) +fi ( 2 ,l) +fi (2, 2 )_ (43) 

In the next three subsections, these results are expressed in terms of the basis integrals. The two-loop fermion pole 
squared masses then follow by plugging eq. (4.3) and the results of section III into eq. (2.51). 



A. Contributions from diagrams with no vector propagators 

The fermion self-energy functions following from the two-loop diagrams of fig. 4 without vector or ghost propagators 



4) 



are: 

vj/ 2,0)J = y ILl y JK: > y KNi y LNj M SF FSF(m 2 ,m 2 K ,m 2 L ,m 2 ,mj 

+V K Nl yL<Njm KK ,m LL M sJTsF (m 2 ,m 2 K ,m 2 L ,m 2 ,m 2 N ) 
+yKNiyL'N' 3 m LL m NN M SF p S p(m 2 ,m 2 K , m 2 L ,m 2 ,m 2 N ) 
+y K ' Nl y LN ' 1 mKK<rn NN ,M SF p SF -(m 2 j , m 2 Ll m 2 Kl m 2 , m 2 N ) 

+\ t:ik yKL>km LL 'M SSF p S (m 2 ,ml m 2 K , m 2 L ,m 2 k ) + X ijk y K ' Lk m K K'M SSF7s (m 2 , m 2 ,m 2 L ,m 2 K , m\) 

+yimy JK3 {\x t]kk Y FS ss(m 2 K , mlm 2 ,m 2 ) + \x ikn A^ fe "V FS sss(m 2 K , m 2 , m 2 im 2 ,m 2 n ) 
+Rc[y LNl y LNj ]V FS sFF(m 2 K , m 2 ,m 2 ,m 2 L , m 2 N ) 
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+Re[y LNt y L ' N ' 3 m LL <m N N ,]V FSsTF (m 2 K , to- , m 2 ,m 2 L ,m 2 N )} 
+yi Kl y JLi y KN3 yLN 3 V S FFFs{ml rn 2 K , m 2 L ,m 2 N , m 2 ) 
+VK'Njy L Nj m KK m LL ,V sTFFS (m 2 ,m 2 K ,m 2 L ,m 2 N ,m 2 ) 
+y KNj y L ' N ' 3 m LL >m NN >V SF -pp S (m 2 , m 2 K , m 2 L , m 2 N , m 2 ) 



+yK'NjyLN'jm to V SF j Ts (m il rn L ,m K ,m Nl m J ) , (4.4) 

iViV'u I 2 2 2 2 2 



^KK' ^NN'\ r 1^2 m 2 _2 m 2 m 2 
LN'j" 

n (2,o)/J = ^V 7 ^ [y^L^m JViv M SFFsF K,TO^m£,TO*,r<) 

+yKNiy L Nj m LLI M SFFSF (m 2 ,m 2 K ,m 2 L , m 2 ,m 2 N ) + y K Ni yLNjm KK ,M SFFSF (m 2 , m 2 L ,m 2 K , m 2 ,m 2 N ) 



+y K Nl y L N 3 m KKI m LL im NN ,M SFTsF -(m 2 ,m 2 K ,m 2 L ,m 2 , m 2 N ) 

+X 3k yKLkM SSFFS (m 2 , m 2 ,m 2 L ,m 2 K , to|) + X ljk y K L k m KK ,m LL ,M ssWFS (m 2 , m 2 , m\ , m 2 K , m 2 k ) 

+y IKi y JK ' j m K K'{\x ijkk Y F sss«^ rnlm 2 , m 2 ) + ^ kn X 3kn V FSSSS (m 2 K ,ml m 2 ,ml,m 2 n ) 

+Rc[y LNt y LNj ]V FSSFF (m 2 K ,m 2 ,m 2 ,m 2 L ,m 2 N ) 

+Rc[y LNt y L ' N ' 3 m LL ,m NNI ]V FSSFT (m 2 Kl to- , m 2 ,m 2 L ,m 2 N )} 



+y IKi y JLi 



NN\r i 2 2 2 2 2\ 

yKNjyLN'jm Vgpppgimi ,m K ,m L ,m N , rrij) 



, L'Nj u / 2 2 2 2 2\ 

+VKNjy J m LL ,V SF p FS {m l ,m K ,m L ,m N ,m j ) 

K'Nj „ 1 2 2 2 2 2\ 

+y J yLN J m KKI V SF - FFS (m l , m L ,m K , m N ,m j ) 

+y K ' Nj y L ' N ' 3 mKK'm LL/ m NNI V SF ^ s (m 2 ,m 2 K ,m 2 L ,m 2 N ,m 2 ) , (4.5) 

where the functions corresponding to each diagram are: 

M SFFSF (x, y, z, u, v) = [(ux -yz- sv)M(x, y, z, u, v) + yU(u, y, x, v) + zU (x, z, u, v) - uU (y, u, z, v) 

-xU(z, x, y,v) - S(x,u, v) + S(y, z, v) + sB(x, z)B(y,u)]/2, (4.6) 
M SFFSF( x >y> z > u > v ) = [{x + u-v-s)M(x,y,z,u,v)-U(y,u,z,v)-U(z,x,y,v) + B(x,z)B(y,u)]/2, (4.7) 

M SFFSF (x,y,z,u,v) = [(x-y-v)M(x,y,z,u,v) + U(x,z,u,v)-U(y,u,z,v) + B(x,z)B(y,u)]/2, (4.8) 

M SFFSF( x >y' z ' u > v ) = [( u ~ V ~ s)M(x,y,z,u, v) + U(x, z,u,v) - U(z,x,y, v)]/2, (4.9) 
M s~ffsf( x ' 2/' z,u,v) = [(x - y - z + u)M(x, y, z, u,v) + U (x, z, u, v) + U(u, y, x, v) - U(y, u, z, v) 

-U(z,x,y,v)]/2, (4.10) 

M sWs: p(x,y,z,u,w) = -M(x,y,z,u,v), (4.11) 



and 



M S sFFs{x,y,z,u,v) = [(v - z - u)M(x,y, z,u,v) + U(z,x,y,v) + U(u,y,x,v) - B(x,z)B(y,u)]/2, (4.12) 

M ssFFs( x 'y' z ' u ' v ) = [( x ~ z ~ s)M(x,y,z,u,v) - U(y,u,z,v) + U(u,y,x,v)]/2, (4.13) 
M SSFFS (x,y,z,u,v) = -M(x,y,z,u,v), (4.14) 



and 



Yf sss(x, y, z, u) = A(u) [A(y) - A(z) + (y - x - s)B(x, y) - (z - x - s)B(x, z)] /2(y - z), (4.15) 

Y FSSS (x,y,y,u) = A{u) [1 + A(y)/y + B(x, y) + (y - x - s)B(x, y')\ /2, (4.16) 

Y FSSS (x,y,z,u) = A(u)[B(x,z)-B(x,y)]/(y-z), (4.17) 

Y FSSS (x,y,y,u) = -A(u)B(x,y'), (4.18) 

and 

^Fssss(x,y,z,u,v) = [(s + x - y)U(x,y,u,v) - (s + x - z)U(x,z,u,v) - I(y,u,v) + I(z,u,v)]/2(y - z), (4.19) 
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^Fssss(x,y,y,u,v) = [{y-x- s)V(x,y,u,v) - U(x,y,u,v) - I(y' ,u,v)]/2, 



FSSSS 
FSSSS 



(x,y,z,u,v) 
(x,y,y,u,v) 



[U(x,y,u,v) - U(x,z,u,v)]/(y - z), 
-V(x,y,u,v), 



(4.20) 
(4.21) 
(4.22) 



and 

VsFFFs{x,y,z,u,v) 
VsFFFs(x,y,y,u,v) 



^ SFFFs( X > f ' Z ' U > V 

^ S FFFs( x ^y'y^ u ' v 

SF~FFs( X > V' Z ' U > V 
V SF F FS {x,y,y,U,V 

^ SFF~F~s( X > y> Z ' U > V 

V 



SFFFS 



(x, y, y, u, v 

^ SF~FFs( X > V' Z > U > V 



(4.23) 



(4.24) 
(4.25) 
(4.26) 
(4.27) 
(4.28) 



~V s -pp FS (x,y,y,u,v) 



SFFFS 



SFFFS 



(x,y,z,u,v) 
(x,y,y,u,v) 



[(s - x + y)(y + u - v)U(x,y,u,v) - yS(x,u,v) + (y + u- v)I(y,u,v) 
+ [A(u) - A(v)][{s -x + y)B(x, y) + A(y)]] /4(y - z) + (y «-» z), 
[(s — x + y)(v — y — u)V(x, y, u, v) + (s — x + 2y + u — v)U(x, y, u, v) 
—S(x, u,v) + (y + u — v)I(y', u, v) + I(y, u, v) 
+ [A(u) - A(v)][(s -x + y)B(x, y') + B(x, y) + l + A(y)/y]] /4, 
[yU(x, y, u,v) - zU(x,z,u,v)]/(y - z), 
U(x,y,u,v) - yV(x,y,u,v), 

[(y + u- v)U(x, y, u, v) + [A{u) - A(v)]B(x, y)] /2(y - z) + [y z), 
[U(x, y,u,v) + (v - y - u)V(x, y, u, v) + [A(u) - A(v)]B(x, y')]/2, 
[(s - x + y)U(x,y,u,v) -(s-x + z)U(x,z,u,v) + I(y,u,v) - I(z,u,v)]/2(y - z), (4.29) 
[U(x, y, u, v) + (x-y- s)V(x, y, u, v) + I(y', u, v)]/2, (4.30) 
{(s - x + y)(y + u- v)U(x,y,u,v) + (u - v + y)I(y,u,v) 
+ [A(u) - A{v)][(s -x + y)B(x,y) + A(y)]}/4y(y - z) + {y <-► z) 
+{2{u — v)[S(x, u, v) + xT(x, u, v)] + u(x + u — v — s)T(u, x, v) 

+v(s + u-v- x)T{v,x,u) + (v- u)[A{x) + A{u) + A(v) -x-u-v + s/4]}/4yz, (4.31) 

y(y + u — v)(x — y — s)V(x, y, u, v) + (sv — su + xu — xv + y 2 )U(x, y, u, v) 

+2{u — v)[S(x, u, v) + xT(x, u, v)] + u(x + u — v — s)T(u, x, v) 
+v(s + u — v — x)T(v, x, u) + y(y + u — v)I(y' , u, v) 
+(v - u)[I(y,u,v) + A(x) + 2A(v) - x - u - v + s/4] 

+ [A(u) - A(v)] [y(s -x + y)B(x, y') + (x - s)B(x, y)+y-u + w]] /4y 2 , (4.32) 

[U(x,y,u,v) - U(x,z,u,v)]/(y - z), (4.33) 
-V(x,y,u,v), (4.34) 



and 



V F ssFF(x,y,z,u,v) 
V FS sFF(x,y,y,u,v) 



V FSS -pp{x,y,y,u,v) 
Vp SSFF (x,y,z,u,v) 
V- 



FSSFF 



(x,y,y,u,v) 
(x,y,z,u,v) 

VFss¥F( x 'y>y' u > v ) 



v 



FSSFF 



Ux — y + s)(y — u — v)U(x,y,u,v) + yS(x,u,v) + (u + v — y)I(y, u, v) 

+[A(u) + A(v)][(s + x- y)B(x, y) - A(y)}]/2(y - z) + {y ~ z), 

Ux — y + s)(u + v — y)V (x, y, u, v) + (x — 2y + u + v + s)U(x, y, u, v) 

+S(x, u, v) + (u + v- y)I{y', u, v) - I(y, u, v) 

+[A(u) + A(v)][(s + x- y)B(x,y r ) - B(x,y) - A(y)/y - l]]/2, 

[(y-x - s)U(x, y, u, v) + I(y, u, vj\ /(y - z) + (y <-> z), 

{x-y + s)V(x, y, u, v) + U(x, y, u, v) + I(y' , u, v), 

[(y-u- v)U(x, y, u, v) + [A(u) + A(v)]B(x, y)] /(y -z) + {y^z), 

(u + v- y)V(x, y, u,v) + U (x, y, u, v) + [A(u) + A(v)]B(x, y'), 

2[U(x,z,u,v) - U(x,y,u,v)]/(y- z), 

2V(x,y,u,v). 



(4.35) 



(4.36) 
(4.37) 
(4.38) 
(4.39) 
(4.40) 
(4.41) 
(4.42) 



Note that for diagrams with the V or Y topology, the limits of identical (or degenerate) squared masses in the second 
and third arguments required separate expressions to avoid the threats of vanishing denominators. Also, the limit 
y — > appropriate for massless fermions is only needed when the corresponding propagator has no mass insertion. For 
the case Vsfff~s> this limit is trivial, since yV(x, y, u, v) vanishes as y — * 0. The remaining non-trivial case involving 



12 



a massless internal fcrmion is 



Vsfffs(x, 0, 0, u, v) = [(s — a;)(w — u)V(x, 0, u, w) + (s — x + u — v)U (x, 0, u, v) — S(x, u, v) 

+1(0, u, v) + 2[A(u) - A(v)} [s - A(x) - xB(0, x)]/(s - x)] /4, (4.43) 

where the function V(x, 0, u, v) is defined in eq. (2.20) of ref. [33] and given in terms of the basis integral functions in 
eqs. (A.11)-(A.13) of that paper. Two useful special cases are: 

Vsfffs(x, 0,0, u,u) = [(s - x)U(x,0, u,u) - S(x,u,u) + 1(0, u,u)]/4, (4.44) 

and for s = x with u ^ v, 

Vsfffs(x,0,0,u,v)\ s=x = {[A(u) -u][A(v) -v] +uv+ (u + v)I(0,u,v)/2}/2(u- v) 

+{2(u -v + x)T(x, u, v) - 3vT(v, u, x) + buT(u, x, v) + 1(0, u, v) 

-3A(u) + A(v) - A(x) +3u-v + 2[A(v) - A(u)]A(x)/x + 3x/4}/8. (4.45) 

The corresponding contribution to the two-loop pole masses [see eqs. (2.51), (2.52), and (4.3)] is given in terms of 
the above results by: 

fi (2,0) = s (2,0)7 + S (2,0)I miI , m Il" /m 2 + mn ,^mi' + n ff m »' + (y IKi y»"y JLJ y"* m * 

+y I ' Ki VJKtV JLj Vi"Ljm u ,m 11 " + y IKl y JKl y J ' L3 y rLj mjj l m II '^B FS (m 2 K , m 2 )B FS (m 2 L ,m 2 )/m 4 I 
+yiK l yjKHy JL] y IL ' 3 m KK 'm LL ,B Ts (m 2 K , m 2 )B FS (m 2 L ,m 2 ) + {2Rc[y I ' Kl yj Kl y JL ^y IL ^ mi rm LLI ] 
+yiK l y JKl [y.jLjyi'L' J m 11 ' m LL ' + y J ' Lj y IL ' j m J j l m LL ^B FS (m 2 K ,m 2 )B FS (m 2 L ,m 2 j )/m 2 I 
+ ([\y IKl \ 2 + |m // ^ / ^|V^]BFs(m^,m l 2 ) + 2Re[m // ,/ K V K 'W^]B 1?s (m^,m2)) 
([| y ^| 2 + \m"" yi ,. Lj \ 2 /m*] [B' FS (m 2 L ,m 2 ) - B FS (m\, m 2 )/m 2 ] 

+2Re[m u »y I " L] y IL ' 1 mLL>} B^ s (m 2 , m 2 )) , (4.46) 

where s — m 2 + ie is taken everywhere on the right side. There is no sum on the index /, but all other indices 
(including I' and I") are summed over as usual. 

B. Contributions from diagrams with one vector propagator 

Next we consider the contributions coming from the two-loop diagrams in fig. 4 that involve exactly one vector 
line. It is convenient to organize these in terms of certain linear combinations of the quadratic Casimir group theory 
invariants for the fermions and scalars appearing in the diagrams, as follows: 

sf' 1)J = yi Kl y JKt gl{llC a (K) + C a (i)-C a (I)}G F s(rn 2 K , m 2 ) 

+ [C a ( K) - C a (i)] G F fs (mj , m 2 K , m 2 ) + C a (I) H FFS (mj , m 2 K , to 2 ) } 

+ {yiK l yj'K' t m JJ 'm KK ' +y I ' Ki y JK ' i m II ,m KK ,} g 2 a {[C a (K) - C a (»)]G 1?Fs (m 2 , m 2 K , m 2 ) 
+C a (I)H FFS (mj,m 2 K ,m 2 )}, (4.47) 

n( 2,l)/J = ^y^ m ^ 5 2|l [Ca(x) + CaW _ Ca(/)]G _ s(m ^ jm2) 

+ [C'a(K) - C a (i)]G F p S (mj,m 2 K ,m 2 )+C a (I) H FFS (mj,m 2 K ,m 2 )} 
(y IKi yj>Kim JJ ' +yp Ki y JKi m II ')g 2 a {[C a (K) - C a (i)]G FFS (m 2 ,m 2 K ,m 2 ) 

+C a (I)H FFS (m 2 ,m 2 K ,m 2 )}. (4.48) 
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The loop integral functions appearing here are: 

G FS (y, z) = [y 2 -(z- sf]M(y, y, z, z, 0) + S(0, y, z)/2 + (3y - z + s)T{y, 0, z)/2 + (y + s)T{z, 0, y) 
+2(y - z + s)T(0, y, z) + sB(y, zf + [3(s -y + z)A(y) + 3(s + y - z)A{z)/2 + 2(y - zf 
-2s(y + z)]B(y, z) + [(s + 3y - z)A(y)/2y + (y + s - z/2)A{z)/z + 6{z-y- s)]B(y, z) 
-33s/8 - 3y/2 + bz/2 + A(y)/2 - A(yf /2y + 2A{z) - 2yA(z)/z + 3A(y)A(z)/z - 3A(zf /2z 
+5msM z -V + s)B(y, z') + (z + y- s)B(y, z)/2 + (1/2 + y/z)A(z) - A(y)/2 + y- s/8] 
+m 2 [(y - z + s)B(y, z') - B(y, z) - A(z)/z - 1], (4.49) 

G Ts (y, z) = 2(y + z- s)M(y, y, z, z, 0) + 2T(y, 0, z) + 2T(z, 0, y) + 4T(0, y, z) + B(y, z) 2 + [2y -2z-2s 
+3(s -y- z)A{y)/y + 3A(z)]B(y, z') + [5A(y)/y + 2A(z)/z - lA]B{y, z) + 3A(y)A(z)/yz 
-3A(y)/y - 2A(z)/z - 6 + 6 m [(s -y- z)B(y, z') + B(y, z) + A(z)/z - 1] + 2m 2 B(y, z% (4.50) 

and 

G F Fs{x,y,z) = [(y - zf - xz + ax - sz]M(0,y, x, z,y) + [(y - zf - (s + x)(y + z)/2]M(0, z, x, y, z) 

-xU(0,x,y, z)/2 + (y- z + s/2 + x/2)U(x,0,y,y) + (z - y)U(x,0, z, z) - yU(y,z,x,y) 

-yU(z, y, x, z) + 55(0, y, z)/4 + S(x, z, z)/2 - S(x, y, y)/2 + (z - y/2)T(z, 0, y) 

+ (s-z + 3y)T{y, 0, z)/A +[(y + z- s)(y - z) - 3{s - z + y)A(z)/4 

+3(s + z - y)A{y)/2]B{y, z') + [(z/4 - y/2)A[z)/z + (s + 3y - z)A{y)/Ay - s/2B(0, x) 

-s/2]B(y, z) + (x + s)B(0, x) + 3A{zf/Az + 3A(y)A(z)/2z - A{yf/Ay + A(x) 

+A(y)/A - (3/2 + y/z)A{z) + 5(y + z)/A - 9s/W 

+ S msK z -V + s)yB(y, z 1 ) + {by -3z + 3s)B(y, z)/2 + 3A(y)/2 + (y/z - 3/2)A(z) +y- s/8}/2 
+m 2 [l + B(y, z ) + (z-y- s)B(y, z') + A(z)/z]/2, (4.51) 

G ffs( x > V' z ) = {4J/ M (°> 2/ ; x > z > y) + {2y + 2z-x- s)M(0, z, x, y,z) + U (0, x, y, z) - 2U {x, 0, z, z) 

-2U(y, z, x, y) - T{z, 0, y) + [B(0, x) + 1 - A(z)/z]B(y, z) - AB(0, x) + I}/ A, (4.52) 

G F p S {x, y,z) = (2x + y - z - s)M(0, y, x, z, y) + (2y - 2z - x - s)M(0, z, x, y, z)/2 + U (x, 0, y, y) 
-U(x, 0, z, z) - U(y, z, x, y) - U(z, y, x, z)/2 + 2T(y, 0, z) - T(z, 0, y)/2 
+[3(s -y- z)A(y)/2y - 3A(z)/2 + y + 3z - s]B(y, z') + [7A(y)/2y - A(z)/2z 
-5/2]B(y, z) + 3A(y)A(z)/2yz - 3A(y)/2y - A(z)/z - 1/2 

+<W( S - y - z ) B (v> z ') + 5B (y> z ) + A ( z )/ Z + l V 2 - m l B iy^ z ')^ ( 4 - 53 ) 

Gpp S {x, y, z) = {2(2s -x + y- z)M(0, y, x, z, y) + (2y - 2z - x - s)M(0, z, x, y, z) + 2U{x, 0, y, y) 
-2U(x, 0, z, z) - 2U{y, z, x, y) - U(z, y, x, z) - 2T(y, 0, z) - T(z, 0, y) 

+ [3 - 2A(y)/y - A(z)/z]B(y, z) + 3}/4, (4.54) 

and 

H FFS (x,y,z) = [(y - zf - xz + sx - sz]M(0,y,x, z,y) + [(x + s)(y + z)/2 - (y - zf]M(0, z,x,y,z) 

+(x + s){z-x- y)V(0, x, y, z) + [(3 + s/x)(z - y) + s}U(0, x, y, z)/2 + (x/2 + y-z + s/2)U (x, 0, y, y) 

+(y- z)U (x, 0, z, z) - yU(y, z, x, y) + ylJ {z, y, x, z) + [sx - 2sy + 3xy + 2y 2 - xz - 2yz]T(y, 0, z)/Ax 

+ {y-z + a)T(0, y, z) + (xy + yz - z 2 + zs)T(z, 0, y)/2x + (3/4 + y/x - z/x)S(0, y, z) - S{x, y, y)/2 

-S(x,z,z)/2+ {x + y- z)I(x',y,z) + (x - y + z)I(x, y, z)/2x + [(y - zf - s(y + z) 

+3(z -y + s)A{y)/2 + 3(y - z + s)A(z)/A}B(y, z') + [(2z - 2y - 3s/2)B(0, x) + {s + 3y - z)A(y)/Ay 

+ (y/2z - l/A)A(z) - 5s/2 - 3y + 3z]B(y, z) + [A(y) - A(z)][2sB(0, x') - (1 + s/x)B(0, x)/2] 

-A(yf/Ay - 3A(zf/Az + 3[A(y) - A(z)]A(x)/2x + 3A(y)A(z)/2z + (z/2x - y/2x + l/A)A{y) 

+ (1 - y/z)(l + z/2x)A(z) - 25s/16 - 3y/4 - sy/8x + y 2 /2x + hz/A + sz/8x - z 2 /2x 

+5 M g[2y(^ - y + s)B(y, z') + (5y - 3z + 3s)B(y, z) + [2y/z - 3)A{z) + 3A{y) + 2y- 3s/ A}/ A 

+m 2 [(y - z + s)B(y,z') - B(y,z) - A(z)/z - l]/2 + (1 - £){( s 2/ - sx - 3x 2 - xy - sz + xz)U(0,x,y, z) 

+2(s + x)(x + y — z)xV(0 1 x, y, z) — 2sx{x + y — z)U'(0, x, y, z) + (s + 2x — y + z)yT(y, 0, z) 

+ {z -y- s)zT(z, 0, y) + (s + y - z)xT(0, y, z) + 2(x - y + z)S(0, y,z) + {y-x- z)I{x, y, z) 
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+2(z -x- y)xl(x',y, z) + 2xs[(y -z + s)B{y, z) + 2A(y) - 2A(z)]B'(0, x) 

+2x(s + y- z)B(0, x)B{y, z) + [A(x) + (s + 3x)B(0, x)\ [A(y) - A(z)\ + (y - z)A{y) 

+(y -2x- z)A(z) + (x - y + z)(y + z - s/4)}/2x, (4.55) 
H F p S (x, y, z) = (2x + y-z- s)M(0, y, x, z, y) + {x/2 + s/2-y + z)M{0, z, x, y, z) + U{x, 0, y, y) + U{x, 0, z, z) 
-4C/(0, x, y, z) - U(y, z, x, y) + U{z, y, x, z)/2 + 8xV(0, x, y, z) + 2T(y, 0, z) + T(z, 0, y)/2 + 2T(0, y, z) 
+ [3{s -y- z)A(y)/2y + 3A(z)/2 -s + y- z}B(y, z') + [7A(y)/2y + A(z)/2z - 25(0, x) - 15/2]B(y, z) 
+3[A{z)/z - l]A{y)/2y - A(z)/z - 7/2 + ^[3 + A{z)/z + 5B(y, z) + (s - y - z)B{y, z')]/2 + m 2 B{y, z') 
+ (1 - £){T(0, y, z) - 2sU'{Q, x, y, z) - 2xV(0, x, y, z) - U(0, x, y, z) + 2[B{Q, x) + sB'(0, x)]B(y, z)}, (4.56) 
Hp FS {x, y, z) = {4yM(0, y, x, z, y) + (x - 2y - 2z + «)M(0, z, x, y, z) - 5U (0, x, y, z) + 2U (x, 0, z, z) - 2U (y, z, x, y) 
+8(x + y- z)V(0, x, y, z) + T(z, 0, y) + [35(0, x) + A(z)/z - l]B(y, z) + 8[A(z) - A(y)]B(0, x') - l}/4 
+6^/2 + (1 - ^{(z -x- y)[xV(0, x, y, z) + sU'(0, x, y, z)] + (z - y)U(0, x, y, z) - yT{y, 0, z) + zT(z, 0, y) 
+ [(s + y- z)B(y, z) - A(y) + A(z)]xB'(0, x) + [A(z) - A(y)][(l + 2x/s)B(0, x) + 2A(x)/s]}/2x, (4.57) 
Hpp S {x, y, z) = {2{2s -x + y- z)M(0, y, x, z, y) + (s + x - 2y + 2z)M (0, z, x, y, z) + 2U (x, 0, y, y) 

+2U(x, 0, z, z) + U(z, y, x, z) - 2U(y, z, x, y) - 4(s + x)V(0, x, y, z) - 2T(y, 0, z) + T(z, 0, y) + 4/(x', y, z) 
+ [1 - 2A(y)/y + A(z)/z]B(y, z) + l}/4 + (1 - + x)V(0, x, y, z) - U(0, x, y, z) - sU'(0, x, y, z) 
-I(x',y,z) + sB(y,z)B'(Q,x)]. (4.58) 

Note that the terms involving functions Gps an d Gps are ^ ne on ^y P ar ts that contribute when the external fermion 
is neutral; they and the functions Gffs, G F -p s , G FFS and G-p-pg are each gauge- invariant and finite in the limit 
s — ► rrij. In contrast, the functions Hffs, B F p S , Hp- FS and Hp-p S are not by themselves gauge- invariant and have 
logarithmic divergences as s — > mf, but they combine with one- loop parts to give a finite, gauge-invariant pole mass. 
This cancellation provides a nice check on the calculations. The resulting contribution to the pole squared mass [see 
eqs. (2.51), (2.52), and (4.3)] is: 

g(2,i) = [| y /*i|2 + \m u 'y IIK tflm*\ g 2 a {\[C a {K) + C a (i) C a (I)]fi(mj, m 2 K , mf) 

+ [C a (K)-C a (i)]/ 2 (mj ,m 2 K ,m 2 ) + C a (I) f 3 (m 2 , m 2 K ,m 2 )} 

+2Rc [y IKl y I,K ' l miFm KK ,] g 2 a {\[C a {K) + C a (i) - C a (I)]f 4 (mj,m 2 K ,m 2 ) 

+[C a (K) - C a (i)}f 5 (m 2 ,m 2 K ,m 2 ) + C a (I)f 6 (m 2 ,m 2 K , m 2 )}, (4.59) 

where: 

fi(x,y,z) = lim G FS (y,z) (4.60) 

S—fX 

= [y 2 - (x - z) 2 ]M(y, y, z, z, 0) + S(0, y, z)/2 + 2(x + y- z)T(0, y, z) + (x + 3y - z)T(y, 0, z)/2 
+(x + y)T(z, 0, y) + B(y, z'){3(x -y + z)A{y) + 3(x + y- z)A{z)/2 + 2[{y - zf - x(y + z)]} 
+xB(y, z) 2 + B(y, z)[(x + 3y - z)A(y)/2y + (2x + 2y- z)A{z)/2z + 6{z-x- y)\ 
+ [1 - A{y)/y]A{y)/2 + [2z - 2y + 3A(y) - 3A(z)/2]A(z)/z - 33x/8 - 3y/2 + bz/2 
+hdv( x -V + Z ) B (V' z') + (y + z- x)B(y, z)/2 + (1/2 + y/z)A(z) - A(y)/2 + y- x/8] 
+m 2 [(x + y- z)B(y, z') - B(y, z) - A(z)/z - 1], (4.61) 
f 2 (x,y,z) = Y\m[G FF s{x,y,z) + 2xGp FS {x,y,z)] (4.62) 

= (x + y- z) 2 M(0, y, x, z, y) + [(y - zf - x 2 }M(0, z, x,y,z) + (x + y- z)U(x, 0, y, y) 

+(z - x - y)U(x,0,z,z) -(x + y)U(y,z,x,y) - yU(z,y,x, z) + 55(0, y, z)/4 - S(x,y,y)/2 

+S(x, z, z)/2 +(x + 3y- z)T(y, 0, z)/A + {2z - x - y)T(z, 0, y)/2 

+ [3{x -y + z)A{y)/2 + 3(z-x- y)A{z)/A + (z-y)(x-y- z)]B(y, z') 

+ [{x + 3y- z)A{y)/Ay + (z - 2x - 2y)A(z)/4z]B(y, z) - A(yf/Ay + 3A(y)A(z)/2z 

+3A{zf/Az + A(x) + A{y)/A - (3/2 + y/z)A{z) + 5(y + z)/A- x/16 

+^[2^(0: -y+ z)B(y, z') + (3a; + 5y - 3z)B(y, z) + (2y/z - 3)A(z) + 3A(y) + 2y - x/A]/A 
+m 2 [(z -x- y)B(y, z') + B(y, z) + A{z)/z + l]/2, (4.63) 
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h(x, y, z) = lim ^H FFS {x, y, z) + 2xHp FS (x, y, z) + 2[B FV (x, 0) + xB Tv (x, 0)]B' FS (y, z) 

+2[B' FV {x, 0) + xB^ v (x, 0)}B FS (y, z) - B FV (x, 0)B FS (y, z)/x} (4.64) 

= (x + y- z) 2 M(0, y, x, z, y) + [x 2 - (y - z) 2 ]M(0, z, x,y,z) + (x + y- z)U {x, 0, y, y) 

+(x + y- z)U(x, 0, z, z)-(x + y)U(y, z, x, y) + yU{z, y, x, z) + [(z -y)/x- 5/4]5(0, y, z) 
-S(x, y, y)/2 - S(x, z, z)/2 + 2(x + y- z)T(0, y, z) + (x 2 + xy - xz - yz + z 2 )T(z, 0, y)/2x 
+ (x 2 — 2y 2 + 2yz — 3xy — xz)T(y, 0, z)/Ax + (3x + y — z)I(x, y, z)/2x + 3(x + y — z)I(x', y, z) 
+ [(5z - 3a; - 5y)A(x)/2x + (x + 3y - z)A(y)/Ay + (2x + 2y- z)A(z)/Az + 6(z-x- y)]B(y, z) 
+ [x(y + z)-(y- zf + 3{y + z - x)A{x) + 3{x -y + z)A(y)/2 + 3(x + y- z)A(z)/A]B(y, z') 
-A(x)A(y)/x - A(y) 2 /4y + (l/x - 3/z)A(x)A(z) + 3A(y)A(z)/2z + (x + 2y- 2z)A(y)/Ax 
-3A(z) 2 /Az + {2xy + 2xz + yz- z 2 )A(z)/2xz + 3x/16 + 3y/8 - y 2 /2x - lbz/8 + z 2 /2x 
+S m [2(y 2 - 3yz + 2z 2 - xy - 2xz)B(y, z') + (5x + 3y - z)B(y, z) - (1 + 2y/z)A{z) + 5A(y) 
+4z - 2y - 3x/ A]/ A + m 2 [(x + y - z)B{y, z 1 ) - B(y, z) - A(z)/z - l]/2, (4.65) 
h{x,y,z) = lim Gp S (y,z) (4.66) 

= 2{y + z - x)M(y, y, z, z, 0) + 2T(y, 0, z) + 2T(z, 0, y) + 4T(0, y, z) + B(y, z) [B(y, z) + 5A(y)/y 
+2A(z)/z - 14] + [2(y - x - z) + 3(x - y - z)A(y)/y + 3A(z)]B(y, z') + 3A(y)A(z)/yz 
-3A(y)/y - 2A(z)/z - 6 + 5 m [(x -y- z)B(y, z') + B{y, z) + A(z)/z - 1] + m 2 2B(y, z% (4.67) 
f 5 (x,y,z) = limJG F p S (x,y,z) + 2G FTs (x,y,z)] (4.68) 

= 2(x + y- z)M(0, y, x, z, y) + 2(y-x- z)M(0, z, x, y, z) + 2U (x, 0, y, y) - 2U (x, 0, z, z) 
-2U(y, z, x, y) - U(z, y, x, z) + T(y, 0, z) - T(z, 0, y) + [l- A{z)/z\ [1 - 3A{y)/2y] 
+ [3{x -y- z)A(y)/2y - 3A(z)/2 -x + y + 3z]B(y, z') + [5A(y)/2y - A(z)/z - l]B(y, z) 
+S m [(x -y- z)B(y, z') + 5B(y, z) + A(z)/z + l]/2 - m 2 B(y, z'), (4.69) 

f 6 {x, y, z) = lim^H F p S (x, y, z) + 2H T p S {x, y, z) + 2[B Fy (a;, 0) + xBp V (x, 0)]B^ s (y, z) 

+2[B' Fy (x, 0) + xB'p V {x, 0)}Bp S (y, z) + B Tv (x, 0)Bp S (y, z)} (4.70) 

= 2(x + y - z)M(0, y, x, z, y) + 2(x - y + z)M(0, z, x, y, z) + 2U {x, 0, y, y) + 2U (x, 0, z, z) 
-2U(y, z, x, y) + U(z, y, x, z) - 45(0, y, z)/x + (1 - 2y/x)T(y, 0, z) + (1 - 2z/x)T(z, 0, y) 
+4T(0, y, z) + 21{x, y, z)/x + 6I(x', y, z) + [3(y -x- z)A{x)/x + 3{x -y- z)A(y)/2y 
+3A(z)/2 + x-y + z]B(y, z') + [5A(y)/2y + A(z)/z - A(x)/x - ll]B(y, z) 
-[1 + A(z)/z]3A(x)/x + 3A{y)A{z)/2yz + (2/x - 3/2y)A{y) + {2/x + l/z)A{z) 
-2(y + z)/x - 1/2 + 5 m [(y -x- 3z)B(y, z') + 3B(y, z) - A(z)/z + l]/2 + m\B{y, z'). (4.71) 

The limits y — > (for massless internal fermions) can be important when there is no corresponding mass insertion: 

/i {x, 0, z) = -{x- z) 2 M(0, 0, z, z, 0) - 5(x - z)U{z, 0, 0, 0)/2 + S(0, 0, z) /2 + xT(z, 0, 0) + xB(0, z) 2 
+ [{x/z - 2)A(z) + 7z/2 - 3x/2]B(0, z) - 3A(z) 2 /z + l\A{z)/2 - l3x/& - 2z 
+5^^ - x)B(0, z)/2 + A(z)/2 - x/8] - m 2 [2A(z)/z + 25(0, z)}, (4.72) 

f 2 {x, 0, z) = (x - z) 2 M(0, 0, x, z, 0) + (z 2 - x 2 )M(0, z, x, 0, z) - xU(0, z, 0, x) + (x - z)[U(x, 0, 0, 0) 
-U{x, 0, z, z) - U(z, 0, 0, 0) /A] + (z- x/2)T{z, 0, 0) - S*(0, 0, x) /2 + 55(0, 0, z)/A 
+S{x, z, z)/2 + [(1 - x/2z)A{z) + x/A - 5z/A]B(0, z) + 3A(z) 2 /2z + A(x) - 13^(z)/4 
+3a;/16 + 2z + S m [3{x - z)B(0, z)/A - 3A{z)/A - x/16] + m 2 [A(z)/z + 5(0, z% (4.73) 

f 3 {x, 0, z) = (x - z) 2 M(0, 0, x, z, 0) + (x 2 - z 2 )M(0, z, x, 0, z) - xU(0, z, 0, x) + (x - z)[U(x, 0, 0, 0) 
+U(x, 0, z, z) - 9U{z, 0, 0, 0)/4] + (x - z + z 2 /x)T{z, 0, 0)/2 - 5(0, 0, x)/2 - S(x, z, z)/2 
+(z/x - 5/4)5(0, 0, z) + (9/2 - z/2x)I{Q, x, z) + [(x/2z - l)A(z) + (3/2 + 5z/2x)A(x) 
-7x/ '4 + 3z/A]B{Q, z) + A(z)[AA(x)/x - 3A(z)/2z - 9/4 - z/2x] - 6A(x) + z 2 /2x - z/8 
+87x/l6 + Sj^Hbx + 3z)B(0, z)/A + 3A{z)/A - 3x/16] - m 2 [A{z)/z + B(0, z)\. (4.74) 
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C. Contributions from diagrams with two or more vector propagators 

We finally turn to the contributions from two-loop diagrams that contain more than one vector (or ghost) propagator. 
Again it is useful to organize the results in terms of common group theory factors. The results for the self-energy 
functions are: 

4 2 ' 2)J = Si&Cail) [gtC b (I)H 1 (m*)+gtC a (G)H 2 (ml)+gll a (K)H 3 (ml,m%)+ , (4.75) 

(2,2)/j = m Uglc a (I) [g*C b (I)H 1 {m 2 I )+glC a {G0H2(m 2 I )+ , (4.76) 

where the required loop integral functions are: 

ffi(x) = -4x 2 M(0,x,x,0,x) + 2S(x,x,x) - (s + 4x)T(x, 0, 0) + 2(s - x)U(x, 0, x, x) + (10x 2 /s - 2s 
-Ax - 4x 3 /s 2 ) ln 2 (l - s/x) + [lls/2 + x - 29x 2 /2s + (4 - 3s/x + Ux/s)A(x)] ln(l - s/x) 
-73s/8 - 5x/2 + (3 + 9s/2x)A(x) + (6/x - s/x 2 )A{x) 2 + S m [(5x 2 /s - s) ln(l - s/x) 
+5x- s/2-sA(x)/x] + (1 -£){4sT(x,0,0) + (3s- llx 2 /s + 8x 3 /s 2 ) ln 2 (l - s/x) 
+ [(2s/x - Ux/s)A(x) - 5s - 6x + 19x 2 /s] ln(l - s/x) + 3s + llx - (s/x + U)A(x) - sA(xf /x 2 
+W0 ~ 5x 2 /s) ln(l - s/x) + sA(x)/x - 5x]} + (1 - £) 2 {sT(x, 0, 0) + (s - x)x 2 /s 2 ln 2 (l - s/x) 
+ [3x -s- 2x 2 /s + (x/s - s/x)A(x)\ ln(l - s/x) + 3s/2 -x + A(x) - sA(x) 2 /x 2 }, (4.77) 

Hi(x) = -2(s + x)M(0, x, x, 0, x) + 4T(x, 0, 0) - 8U(x, 0, x, x) + 42 - 32A(x)/x + 10A(x) 2 /x 2 
+ (14 - 6x/s)(l - x/s) ln 2 (l - s/x) + [24(1 /x - 2/s)A(x) - 36 + 52x/s] ln(l - s/x) 
+5 m 8[A(x)/x - 1 + (1 - 2x/s) ln(l - s/x)] + (1 - £){-4T(x, 0, 0) + 8A(x)/x 
-2A(xf/x 2 - 6 - 2x/s + (8x/s - 6 - 2x 3 /s 3 ) ln 2 (l - s/x) + [(14/ s - 8/x)A(x) + 12 
-12x/s - 4x 2 /s 2 ] ln(l - s/x) + 6 m 2[l - A(x)/x + (2x/s - 1) ln(l - s/x)}} 
+(1 - £) 2 {-T(x, 0, 0) + A(xf/x 2 - A(x)/x + x/s - 3/2 + [(1/x - l/s)A(x) - 2x/s 
+2x 2 /s 2 } ln(l - s/x) + (x 3 /s 3 - x/s) ln 2 (l - s/x)}, (4.78) 

H 2 (x) = (x - s)(s + 2x)M(0,0,x,x,0) + 2x 2 M(0, x, x, 0, x) + (x - s)U(x,0,x,x) - (s + 2x)T(x,0,0) 
-S(x,x,x) - A(x) 2 /x+ (7 + 6s/x)A(x) - 16x - lls/2 + (1 + 2x/s)(x - s)ln 2 (l - s/x) 
+(x - s)[(3/s + l/x)A(x) - 12x/s - 7] ln(l - s/x) + <%jgl3s/4 + (1 - £){s(s - x)M (0, 0, x, x, 0) 
+sT(x, 0, 0) + 2sA(x) 2 /x 2 - 5(1 + s/x)A(x) + 5x + lis + (1 - x/s)[(4x + 3s) ln(l - s/x) 
+5(1 + s/x)A(x) - 5x - 6s] ln(l - s/x)}/2 + (1 - £) 2 {-sT(x, 0, 0) + (1 + s/2x)A(x) - 9s/4 - x/2 
+ (1 - x/s)[(3s + x)/2 - (x + s) ln(l - s/x) - (1 + s/x)A(x)] ln(l - s/x)}/2, (4.79) 

H 2 (x) = 3(s - x)M(0, 0, x, x, 0) + (s + x)M(0, x, x, 0, x) + 7T(x, 0, 0) + 4U(x, 0, x, x) + A(x) 2 /x 2 
-64A(x)/3x + 112/3 + (1 - x/s)[81n(l - s/x) + 9A(x)/x - 85/3] ln(l - s/x) 
+5 m [(x/s- l)ln(l - s/x) - A(x)/x - 7/2] + (1 - £){(x - s)M(0,0,x,x,0) - T(x,0,0) 
-2A(x) 2 /x 2 + 10A(x)/x - 16 + (1 - x/s)[(x/s - 3) ln(l - s/x) + 11 - 5A(x)/x] ln(l - s/x)}/2 
+ (1 - €f{T(x, 0, 0) - 3A(x)/2x + 11/4 + (1 - x/s)[ln(l - s/x) + A(x)/x - 5/2] ln(l - s/x)}/2, (4.80) 
H 3 (x, y) = { [-4s(s - x) 2 + 40s 2 y + 4sxy + 4x 2 y + 48(x - 2s)y 2 ]yT(y, y, x) + 2x(s - x) [9y(s + x) 

-(s - xf - 36y 2 ]T(x, y, y) + [-2(s - x) 3 + 22s 2 y - 8sxy - 14x 2 y + 48(x - 2s)y 2 ]S'(x, y, y) 
+ (s - x) 4 B(0, x) + [2s 2 + 2sx - 4x 2 + 24(x - 2s)y] A(y f + [4(s - x) 2 - 24(s + x)y]A(x)A(y) 
+ [112sxy - 4x(s - x) 2 - 92s 2 y + 28x 2 y + 96(2s - x)y 2 ]A(y) + [(s - x) 3 - 4s 2 y - lOsxy 
+ Ux 2 y + 24(s - 2x)y 2 ]yl(x) + (5s/ 4 - 2x)(s - x) 3 - 35s 3 y/2 + lls 2 xy/2 - llsx 2 y - 10x 3 y 
+8s 2 y 2 + 50sxy 2 + 20x 2 y 2 + 96(x - 2s)y 3 }/15y(s - x) 2 + Sj^Ss/4, (4.81) 
H 3 (x, y) = { [6(s - x) 2 + 12xy + 16y 2 - 28sy]T(y, y, x) + 20x(x - s)T(x, y, y) + (20x - 20s + 16y)S(x, y, y) 
+ (8 + 6x/y - 6s/y)A(y) 2 + WA(x)A(y) + (10s - lOx + 8y)A(x) + (48s - 64x - 32y)A(y) 
+ 19s 2 /2 - 53sx/2 + 17x 2 - 34sy + 8xy + 32y 2 }/3(s - x) 2 + 2(1 + A(y)/y)B(0, x) - 5 m , (4.82) 
H 4 (x,y) = {[48(2s - x)y 2 - s(s - x) 2 - 20s 2 y + 16sxy - 44x 2 y]yT(y,y,x) + x(x - s)[(s - x) 2 /2 

-12y(s + x + 6y)]T(x, y, y) + [y(13s 2 - 2sx - llx 2 + 96sy - 48xy) - (s - x) 3 /2]S*(x, y, y) 
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+ [23s 2 - 22sx - x 2 + 48sy - 24xy]A(yf + [(s - xf + 24(s + x)y]A(x)A(y) + [(s - xf/A 
+(x - s)y(s + llx) + 2A(2x - s)y 2 ]A(x) + [-x(s - x) 2 + 96{x - 2s)y 2 + 22s 2 y - 92sxy 
+22x 2 y]A(y) + (5s/16 - x/2)(s - xf + 96(2s - x)y 3 - 70xy 2 (s + x) - 25s 3 y/4 + 23s 2 xy/2 
+l9sx 2 y/4 - I0x 3 y + 62s 2 y 2 }/15y(s - xf + (s - x) 2 B(0, x)/60y, 
H 4 (x, y) = { [3(s - xf - 8sy + 24xy - l6y 2 ]T(y, y, x) + I6x(x - s)T(x, y,y) + 16(x-y- s)S(x, y, y) 
+ (12x/y - 12a /y - 8)A{yf - 16A(x)A{y) + 8{s-x- y)A{x) + 8(3s - x + 4y)A(y) + Is 2 
-20sx + I3x 2 - 38sy + Mxy - 32y 2 }/3(s - xf + [1 + A(y)/y]B(0, x). 

The corresponding contributions to the pole squared mass are therefore: 

fif' 2) = g 2 a C a (I) [gtCimim 2 ) + g 2 a C a (G)F 2 (m 2 ) + g 2 J a (K)F 3 (m 2 ,m 2 K ) + g 2 a I a ( t )F 4 (m 2 ,m 2 )} 

where the required functions can be written compactly in terms of logarithms and dilogarithms: 

Fi(a;) = lim{2ffi(a:) + 2xH 1 (x) +4[B Fy (a;,0) + xBp V (x,0)}[B' FV (x,0) + xB^ v (x,0)} 
-[B FV (x,0)} 2 /x + x[Bp V (x, 0)] 2 } 
= x Ul/4 + IOtt 2 - 27Ea; + IShfx + 24((3) - 16tt 2 In 2 + <%jgl2(ma; - 1)1 , 
F 2 (x) = lim\2H 2 (x) + 2xH 2 (x)] 

= x [l093/12 - 87r73 - (179/3)mx + llhi 2 a; - 12£(3) + 8tt 2 In 2 + <%jg(3/2 - 2lnx) 
F 3 (x,y) = lim [2H 3 {x) + 2xH 3 {x)} 

s—>x 

= -37.T/3 - 12y + (26x/3)hix - 2xhfx + 4yln(x/y) - 2{y 2 /x) \n 2 (x/y) 
+8(x + y)f(VyJ^) ~ 4(z + y 2 /x)[U 2 {\ - v/x) + tt 2 /6] + 6 m x/2, 
F 4 (x,y) = lim \2H 4 (x) + 2xH 4 (x)] 

= -125x/12 + Uy + (19a;/3)Lna; - xhfx - 6y]n(x/y) + (y 2 /x) \n 2 (x/y) 
+8(x - y)f(VyJ^) + Hv 2 /* - x)[Li 2 (l - y/x) + tt 2 /6], 

I 



(4.83) 



(4.84) 



(4.85) 



(4.86) 
(4.87) 
(4.88) 

(4.89) 
(4.90) 

(4.91) 
(4.92) 

(4.93) 



with 



f(r) = r{Li 2 ([l-r]/[l + r]) 

-Li 2 ([r- l]/[l + r]) +tt 2 /4}. (4.94) 



Here are some useful special cases: 
'26 T 



F 3 (x,0) 



F 3 (x,x) 



F 4 (x,x) 

F 4 (x,0) 
F 3 (0,y) 



x 



3 



-\nx — 21n x 



37 

y 



47T 2 

"3" 



+?)■ 

f 26— —2 73 8tt 2 

x[ — mx — 21n x 1 

V 3 3 3 



MS 

2 

/43 



19, 



Hii + f ln *- lA )' 

19— —2 125 

_ ma; _ m . T ___ 

F 4 (0,y)=0. 



27T 2 

3 



(4.95) 

(4.96) 
(4.97) 

),(4.98) 
(4.99) 



In the MS scheme with = 1, the expressions of 
eqs. (4.87) above agree with those found originally in [16] 
(see also [18]). In the DR scheme with = 0, the same 
equations are in agreement with the results of ref. [17]. 
In particular, the function F 3 (x,y)/2x found here was 
given in three different mass expansions in eqs. (18)-(20) 
of [17]. 



V. EXAMPLES AND APPLICATIONS 

In this section, I present some applications of the pre- 
ceding general results. These are all taken from the min- 
imal supersymmetric standard model (MSSM), and in- 
clude all effects at one-loop order, but only terms that in- 
volve the strong gauge coupling constant in the two-loop 
part. All of the couplings and masses appearing below 
are tree-level running DR parameters in the MSSM with 
no particles decoupled. (This means that contributions 
listed above containing <5j^g and m 2 are not present.) The 
conventions used here for these couplings and masses are 
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identical to those found in section II of ref. [56] and sec- generations is neglected. Throughout the following, each 

tion II of ref. [57], and will not be repeated here for the index that appears on the right-hand side of an equa- 

sake of brevity. Note that the formalism allows for ar- tion but not on the left-hand side is implicitly summed 

bitrary CP violation within the MSSM, and takes into over. The name of a particle is used in place of its renor- 

account sfermion mixing within each generation, but for malized, tree-level squared mass when appearing as the 

simplicity the effects of possible sfermion mixing between argument of a loop integral function. 

I 



A. SUSYQCD corrections to the gluino mass 

In this section, I will study the specialization of the above results to the case of the gluino mass. (The results below 
partly overlap with recent independent results of Y. Yamada in [32].) Applying the formulas of sections III and IV, I 
find that the two-loop gluino pole mass, including all SUSYQCD effects, can be written as: 

M 2 - iT g M g = m 2 g + Ifif + ^fif) , (5.1) 

where m g is the tree-level running gluino mass (often seen as M 3 in the literature), and 

n^ 2) = SI q (2C q - C G )^Lq j R\.L\ k Rq k M S FFSF{qj-,q,qAk,g) -Re[L qk Rl k ]m q m g M SF p SF {qj,q,q,q k ,g) 

+ \{\ L q 3 L qA 2 + ^Rq^WghigFFg-piq^q.qAk.g) + | Lfc Rq J 2 VTi\ Mgjrpgp % , q, q, q k , g) 

-BelL^R^mqmgMgppgpiqj^, q, q k ,g) + Re[Lq.i?|. J L^i?|Jm^m?M sFFsF (^, q, q, q k ,g)} 

+^ I qCq{(\ L q j L q k \ 2 + I Rq, R<j k 1 2 ) ^sfffs % , 9, 9, 9, 5fe ) - 2Re[L qj R\^m g m- g Vgpppg (qj , q, q, g, q k ) 
+2Re[Lg j Rl.Ll k Rq k ]m 2 g V SF pp S (q j ,q,q,g,q k ) + (\L qj R qk \ 2 + |i?g 3 iqJ 2 )m2v sWFS (^,g,g,5,q fc ) 
-2Re[Lg k R* Sk ]m g m g V SF pp S (qj, q, q, g, q k ) + 2Re[L Sj Rt j L 9k RtJm^mp^ SFTFs (q j ,q, q, g, q k ) 
+Vfssff {q, qj , q 3 ,q,g) + 2Re[L gj R*~ k (R*~, L gk + L*-. R qk )}rn 2 q m 2 g Yp S gpp(q, q ,q k ,q, g) 
-2Re[L qj Rl.]m q mgWp SSF p{q, qj,q jy q, g) - 2Re[L gj R\ j ]m g m g W FSSF -p{q, q 0l q k , q,g) 

+ \ I L v L t - R q 3 R t 1 2y FSSS (q, q~j , qj , qk ) 

-Re[L*-.R gk {L gj L*- n - R gj R*- n ) (L*- k L qn - R* gk Rq n )]mgm g Y 7sss (q, q h q k , q n )| 

+4I 2 [4{B FS (q,qj) - 2Re[L 9j R^ j ]m g m s B Fs {q,q j )}{B' FS (Q,Q k ) - 2Re[ J L^ J R^]m Q m § B^(Q, Q k )} 
-B F s(q,qj)^Fs(Q,Qk)/m 2 g + 4m g m Q Re[L qj Rl j }Re[LQ k R^ k }Bj: s (q,q j )Bp S (Q,Q k ) 

+ {AC q - 2C G )I q {f 1 (g,q,q j ) - 2Re[L qj R* qj ]m q mgf 4 (g,q,qj)} 
+AC G I q {f 3 (g, q, qj) - 2Rc[L qj R q [ j ]m q m g f 6 (g, q, q )} 

+C 2 G [Fx (g) + F 2 {g) + F 3 (~g, g)] + C G I q [2F 3 (g, q) + F 4 (g, &)] . (5.3) 

Here C q = 4/3, C G = 3, and I q = 1/2. The symbol q is summed over the 6 symbols (u,d,c,s,t,b), and the squark 
squared-mass eigenstate labels j, k, n are summed over 1, 2. In terms where Q appears, it is also independently summed 
over (u, d, c, s, t, b). The symbols L qj and R qj describe the squark mixing and CP violation; they denote the left-handed 
and right-handed squark content amplitudes of each squark mass eigenstate, as defined in ref. [56]. The loop integral 
functions were listed in sections III and IV above, in terms of basis functions that can be evaluated numerically using 
[43]. Each of them should be evaluated with s = m~ (the tree- level squared mass), with an infinitesimal positive 
imaginary part. A computer program implementing this result is available from the author on request. 

As a non-trivial check, I have verified that this result for the gluino pole mass is invariant under changes in the 
renormalization scale governed by the two-loop SUSYQCD renormalization group equation for the running gluino 
mass [58, 59], [39], up to consistently neglected terms at three- loop order in SUSYQCD and two- loop order in the 
other couplings. 
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In the limit that squark mixing and quark masses are neglected, the expressions above simplify and can be given 
analytically in terms of polylogarithms [60] . The result for the pole squared mass is then: 

M?-iT- g M- g = m l + ^[c G ^ a \g) + i q nf h \9,q j )] 



+ . 5 3 



(16tt 2 ) 2 



C 2 G nf a \g) + C G I q uf b \g, qj) + C q I q uf' c \g, q ) + I 2 q nf d \g, q 3 ,Q k )\ , (5.4) 



where 



Uf' a) (x) = ar[10 — 6biar], (5.5) 

Uf' b) (x,y) = 2x[(l-y/x) 2 ln(l-x/y)+T^y-2 + y/x], (5.6) 

Ilf' a \x) = a;[77 + 107r 2 + 12C(3)-87r 2 ln2-78l^a; + 27hi 2 a;], (5.7) 

Uf' b) (x, y) = (2Ay - 18a: - ±y 2 /x)U 2 (l - y/x) + Ax(x - y)M(0, y, y, 0, x) + 4(x 2 - y 2 )M(0, x, y, 0, y) 

+2(x - y) 2 [M(0, 0, y, y, 0) + 2M(0, 0, x, y, 0) - (A/x) ln 2 (l - x/y)] + 8(y - x)f(^y~Jx~) 
+ [{12y 2 /x - 12y - 2x)]^x + (2Ay - 16x - 6y 2 /x)]ny + Alx - AOy - y 2 /x] ln(l - x/y) 
+ (10y -8x- 2y 2 /x) ln 2 (x/y) - 18a;lna;lm/ + (41a; - Uy)kiy + (9a; + 20y)rna: 
+2y(l - y/x)n 2 /3 - Alx - 9y, (5.8) 
ffi 2 ' c) {x, y) = 8x(y - x)M(0, y, y, 0, a;) - 4(x - y) 2 M(0, 0, y, y, 0) + (4a; - 24y + 24y 2 /x)Li 2 (1 - y/x) 



-2{x - y) 2 (x/y 2 - 2/x) ln 2 (l - x/y) + [(12a; - 8y + 4x 3 /y 2 - 8x 2 /y)hix 

+ {Ax -8y- Ax 3 /y 2 + 8x 2 /y)h\y + 4(y - x)(x/y + 7y/x)] ln(l - x/y) 

+ [12y 2 /x - 12y + Ax 2 /y - 2x 3 /y 2 + 2x] In 2 (x/y) + [Ax 2 /y - 2Qx + 2Ay) \n(x/y) 

+ (Ax 2 /y - 2x 3 /y 2 - Ax/3 - 8y/3 + Ay 2 /x)n 2 - 22x + Ay, (5.9) 



Iif d \x,y,z) = n^''' ) (a;,y)[(3 + 5z/a;)(l-z/a;)ln(l-a;/z) + 31nz-5z/a;-2]/2. (5.10) 

The integral M(0, x, y, 0, y)\ s =x can be reduced using recurrence relations to results found in refs. [18, 19], and 
was given in the present notation in [43]. The analytic formulas for the master integral cases M(0, 0, x, y, 0) and 
M(0, y, y, 0, x)\ s=x were also given in terms of polylogarithms in [43]. The integral M(0, 0, y, y, 0) was originally found 
in [48], and listed in the present notation in [42]. The function / appearing in eq. (5.8) was defined in eq. (4.94) 
above. 

To illustrate the numerical significance of the two-loop correction, in fig. 6 I show the fractional difference between 
the real part of the gluino pole mass and the running renormalized mass (evaluated at a renormalization scale equal 
to itself) as a function of the ratio of the squark masses (assumed degenerate) to the tree-level gluino mass. Also for 
simplicity, the top quark mass is neglected. In most realistic models of supersymmetry breaking, most of the squark 
masses are larger than about 0.8mg. Then the two-loop contribution to the gluino pole mass is positive, and from 1% 
to 2% for comparable gluino and squark masses, but it is larger when m q ^> m§. 

The scale-dependence of the calculated pole mass is shown in fig. 7. To make this graph, a reference renormalization 
scale Qo is chosen such that the running gluino mass evaluated there is equal to it, i.e. Qo = m g (Qo). Then, for three 
different values of the ratio mq/m g at the scale Qo, namely 0.9, 1.5 and 3.0, the one-loop and two-loop gluino pole 
masses are computed as a function of Q. To do this, the relevant running parameters m g , m q , and as = g 2 /A-K are 
evolved using their two-loop renormalization group equations from Qo to Q, and then the pole mass is recomputed 
using eqs. (5.4)-(5.10). In the ideal case, the lines shown would be exactly horizontal. The scale dependence of the 
one-loop and two-loop results is about the same if the gluino mass is less than or about equal to the squark masses (as 
for example in gaugino-mass dominated or "no-scale" models). For heavier squark masses, it is significantly improved 
by going to two-loop order. Note that as usual, the scale-dependence of the one-loop approximation is considerably 
less than the difference between the two-loop and one-loop results. This strongly suggests that the scale-dependence 
should not be used as an indicator of the accuracy of the two-loop approximation either. A naive estimate of the 
size of the three-loop SUSYQCD contribution to the gluino pole mass can be obtained by considering the cube of the 
one-loop fractional contribution, and so is perhaps of the order of a few tenths of a per cent. 

As mentioned above, results equivalent to those above in the limit of no squark mixing have previously been obtained 
independently by Y. Yamada in [32] , where numerical results were given and the analytical form for the pole mass of 
the gluino was shown explicitly in the limit rriq 3> trig. I have checked that the results found here do agree with those 
in ref. [32]. (When comparing the numerical results, it is useful to note that small differences, formally of three-loop 
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FIG. 6: The SUSYQCD corrections to the gluino pole mass M§, in the simplifying approximation of degenerate squarks, 
with no squark mixing or quark masses and Yukawa and electroweak effects neglected for simplicity. Here m squar k = rriq 
and m g i u ino = mg (often seen in the literature as M3) are the tree-level running DR' squark and gluino mass parameters, all 
evaluated at a renormalization scale Q = rrig(Q), and M g i u j nopo i c = M§ is the square root of the real part of the pole mass. In 
the first panel, the dashed line is the one-loop result, and the solid line is the two- loop result, for the ratio of the gluino pole 
mass to the running gluino mass evaluated at itself. The second panel shows the difference between the two-loop and one-loop 
results. The computations were done by specializing eqs. (5.4)-(5.10) in the text, using as{Q) = gf/4-7r = 0.095. 

order, arise due to the fact that the present paper works in terms of perturbative corrections to the pole squared 
mass, while ref. [32] computed results for the pole mass.) 

Further accuracy can be obtained by including the contributions of Feynman diagrams that are of order as times 
Yukawa or electroweak couplings squared. Such corrections could be particularly important if the gluino is relatively 
light, since we are working in a non-decoupling scheme. These results are implicitly contained above in section 
IV; obtaining their explicit form is only a matter of plugging in the couplings and masses of the MSSM. I will not 
do that here, because this paper already contains plenty of lengthy formulas, but only note that it can be done 
straightforwardly by use of a symbolic manipulation program, for example. 



B. SUSYQCD corrections to quark masses in minimal supersymmetry 

As another application of the general results above, consider the relation between the running and pole masses of 
the quarks in the MSSM, found earlier in [29-31]. Using the results of the present paper, I find for the top quark: 

Mi - iT t M t = mf + J^fiN + jS^^ (5 - H) 

where: 

g« = g 2 3 C q {m 2 (10 - 6km?) + 2B FS (5, t ) - 4Re[L f . J R|. ]m t m§B ?s (g, i 3 )} + 2|Y^o \ 2 B F s (t, $) 

+2Rc[(y tT0 o) 2 ]B 1?s (i, + (|F Ifc0 + 1 2 + \Y- bt4> - \ 2 )B FS {b, 0+) + 2Tte[Y JH +Y lt(f> -)rn b mtB- Ps (b, <f>+) 
+(|1Wj I 2 + \Y-tNdi \ 2 ) B fs(N u tj) + 2Re[Y t ^Y mjj }m t m^B Ts (N u t 3 ) 
+ (\Y t cr b * | 2 + \Y Wiij \ 2 )B FS (C^) + 2Re[Y teii ,Y ieilj }m t m d B Ts (C i , b 3 ) 

+ ie 2 m 2 [10 - 6lKm 2 ] + [(g 2 + g' 2 )/i - 2Q t T*g' 2 + 2Q 2 t g' A /{g 2 + g l2 )]B FV {t, Z) 

+2Q t9 ' 2 [Q t g' 2 /(g 2 + g' 2 ) - T 3 *]B w (t, Z) + g 2 B FV (b, W)/2, (5.12) 
h (2) = C q (8C q -AC G ){L i .Ri,LlR ik M SFFSF (t j ;~gr9A,t) + 
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FIG. 7: The scale-dependence of the computed gluino pole mass M g i u i nopo i c = M§, in the simplifying approximation of 
degenerate squarks with no squark mixing or quark masses and Yukawa and electroweak effects. Here the top pair of lines are 
for rriq(Qo) /mtj(Qo) ~ 3, the middle pair for m,q(Qo) /m g (Qo) = 1.5, and the bottom pair for rriq(Qo) /m g (Qo) — 0.9, where 
m g and m,- are the running gluino and squark mass parameters evaluated at a reference renormalization scale Qo = m g (Qo), 
with as(Qo) = <?I/47r = 0.095. In each case, the solid line is the two-loop result, and the dashed line is the one-loop result. 
The parameters m g , rriq, and as are each run from the reference scale Qo to the scale Q using two- loop renormalization group 
equations, and the gluino pole mass is then recomputed. The computations were done by specializing eqs. (5.4)- (5. 10) in the 
text. 
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2 + I^^J 2 )^M s ^ SF (^,g,5,t fe ,t) 
^{Ll j ^ J ]mtrrigM SF - ¥s - w {i j ,g, g, i k ,t) + Re[L fj i?|.L ffc i2|Jm?m|M sl .p sF (^,p, g, i k ,t)j 



+8C q I q ^V S FFFs(tj,g, 9, q, Qk) + 2Re[L i M~ j L* 9k Rg k ]m t m q V SFFTs (i : j,g, g, q, q k ) 



1 



- 2R 4 L t j R i j ] m t™>gV S FFFs( t j> 9, 9, 9, Qk) + 2 m l V SFFFs(*J > 9, 9, 9, Qk) 
-2Re[L Sh RtJm q m,gV SFFFS (tj,g, g, q, q k ) + 2Re[L ij R^L qk R q Jm q m t mp/ SFFFS (t j ,g, g, q, q k )j 
-^Cl{yFSSFFi.g,tj,t ,t,g) - 2Re{L i] R^]m t m g [V FSSFF (g ,tj ,tj,t,g) + V FSSFF (g,tj ,tj,t,g)} 



+2(Re[L~ t] R* L ik R* ik ] + |L f , J^/Km^^s, t 3 , i k , t, g)} + ^.L} - Ri^fYpsssM^h) 



-Re[LlR ik {L h Ll - R~ t] B& (L L L* h - R in RIJ]m t m- g Y FSSS (g,tj,i k ,t n ) 
+{B FS (g 1 t J ) - 2Re[L {i Ri.]mm g B FS (g,t j )}{B' FS (g,t k ) - 2Re[L fjk i^Jm t m g B^ g (ff,i k )} 
+Re[L- t .Ri j }Re[L h Rl}mp FS (g, tj)B FS (g, t k ) - \L~ tj R ik \ 2 B FS {g, ij)B FS (g, t k )/m 2 } 

+C G C q h{t, g, ii) + 2{C G - C q )C q f 2 (t, g, ij) + 2C 2 f 3 (t, g, ij) 
~2C q Re[L ij Rl}m t m g {C G .U(t, g, ij) + 2{C G - C q )f 5 (t, g, ij) + 2C q f 6 (t, g, ij)} 
+C%F 1 (t) + C G C q [F 2 (t) + F 3 (t, g)] + C q I q [2F 3 (t, q) + F 4 (t, q 3 )] . 



(5.13) 



The gauge group constants are C q = 4/3, C G = 3, I q = 1/2, Q t — 2/3 and Tg = 1/2. The results for the bottom 
quark can be obtained by taking t <-» b everywhere, with Q t — » Qb — —1/3 and Tg — ► Tg = —1/2. The one-loop part 
given by eq. (5.12) was given in a different notation in [61]. It includes all effects, including those due to the exchange 
of virtual neutral Higgs scalars [</>° = (h°, H°, G°, A ) for j = 1,2,3,4] charged Higgs scalars [<pj = (G ± ,i? ± ) for 

j = 1,2], neutralinos [Nj for j — 1, 2, 3, 4], and charginos [Cj for j = 1, 2] and top and bottom squarks. In each of the 
loop integral functions in eqs. (5.12) and (5.13), one should take s = m 2 with an infinitesimal positive imaginary part. 
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[Here I have only listed the pure SUSYQCD corrections explicitly in the two- loop part, but the corrections involving 
Yukawa couplings and scalar trilincar couplings are also given implicitly by specializing the results of section IV.] 

I have checked that the result given above is consistent with the two-loop renormalization group equations, in the 
sense that the pole mass is invariant under changes in the renormalization scale given by the two-loop SUSYQCD 
renormalization group equation for the top quark mass. As another non-trivial consistency check, I have verified 
that in the (clearly unrealistic) supersymmetric limit, the top-quark pole mass given above is precisely equal to the 
top-squark pole mass as found in eqs. (5.28)-(5.30) of ref. [34]. [In the published and original preprint versions of that 
paper, the term in eq. (5.30) proportional to I q was missing a factor of 2.] 

The two- loop SUSYQCD contribution and Yukawa contributions had already been found in refs. [29-31], using a 
method in which loop integrals are evaluated using an expansion in small mass hierarchies. In principle, the present 
paper generalizes this somewhat, since here I use two-loop integral basis functions without mass expansions. However, 
as was recently pointed out in [31], the actual top and bottom quark masses are such that the leading terms in the 
mass expansion already give an extremely good approximation throughout most of the parameter space left available 
to supersymmetry. I have also checked agreement with eqs. (57)-(62) in ref. [29]. References [29-31] also include an 
extensive study of the impact of the two-loop top and bottom quark mass corrections in the MSSM. 

C. SUSYQCD corrections to neutralino and chargino masses 

In this section, I present the two-loop corrections to neutralino and chargino masses that involve the strong coupling. 
These arise from gluon and gluino propagator lines added to the one-loop Feynman diagrams that involve quarks and 
squarks, and so are parametrically of order asy 2 , asy 2 , asg 2 , asgg' and asg' 2 - These two- loop contributions are 
evaluated by specializing the results above, and do not require the neglect of W and Z boson masses. 

Using the general results of sections III and IV, one can write the neutralino pole masses as: 

Ml - iTjj. Mjj. = m% + ^Hg + p^ng, (5.14) 

for i = 1,2, 3, 4, where m 2 ^ are the tree-level squared mass eigenvalues, and the one-loop part is [61, 62], 

S S = 2n f(\ Y fNj; I' + \ Y 7nJ 3 \ 2 ) b fs(L Ij) + ^n f Re[Y f ^ f -,Y mfj ]m^m f B Ts (f, f 3 ) 
+2|%^. o| 2 B FS (A> J ,0O) + 2Rc[(Y^. „) 2 ]m^ m^B FS {N 3 ^l) 

+ 2 (\ Y c+Ni4>-\ 2 + \ycjN^ k \ 2 )VFs{C 3 A + k ) +4Re[Y e +^ -Y 6 -^ +]m^m^.B 1?s (C' J ,^+) 
+2(g 2 +g' 2 ) {|0^| 2 B Fy (A> ,Z) - Rc[(0^) 2 ]m^ m^B w (iV,, Z)} 

+2g 2 (\0^\ 2 + \0*\ 2 )B FV (C 3 ,W) + 4g 2 R C [O^Of;}m^m dj B Tv (C 3 ,W), (5.15) 
and the two-loop SUSYQCD part (i.e., the part involving the strong interactions) is 

n f t = 16 -93 2 {^l^| 2 + \Y- q N^\ 2 )h(^,q,q 3 ) + Re[Y q ^.,Y q ^]m^ m q f 4 (N t ,q,q 3 ) 

-2Re[Y qfl .-. Y'^R^LqjMsFFSFiqj, q, Q, <?fc,s) - 2Re[Y ?J y. -» Y*^^ L qj Rq k ]m 2 M s -pp SF (qj , q, q, q k , g) 

+ R 4Y q m q * Y q N t q' k L V L h + Y qN^ Y qmq k R l R V m Ni m 9 M SFFSF(qj > 9, 9, Qk, ff) 

-2Re[Y qfr ..,Y q ^,-,L qj R qk + Y-^Y-^ k RlL* q Jm^m q M SFTsF {q 3 , q, q, q k , g) 
+2Re[Y q „ iq , Y* At . h R qj R qk + Y*^ Y q ^., L qj L qk }m q m g M SFFSF (q 3 , q, q, q k , g) 
+ R 4Y q N ir .Y q ^~,R qj R qk + Y-ft,~. Y-^.- k V q . L-Jm^. mgm 2 M SFFSF (q 3 , q, q, q k ,g) 

+ (\ Y q m q * I 2 + \ Y qN iqj \ 2 ) Y FSSFF(q, qj,q 3 ,g, q) + 2Re[Y^ Y^~Jm^m 9 V ?SSFF (g, Qj, 9j,9, q) 
-2Re[(Y q ^..,Y*^ q , + Y ^Y-^ qk )R qj L qk }m q m g V FSSFF (q, q 3 ,q k ,g, q) 

-2Re[Y q ^.-, Y-fj.- k (R qj L\ k + L qj R\ k )]m# t m 2 q m- g V FSSFF {q, q 3 , q k , g, q) 

+(\ Y q N iq * L Qk I 2 + \ Y q N iqj R qA 2 )VsFFFs{q 3 ,q,q,gAk) + 2Re[Y g ^.-,Y^.-Jm^.m g V SFFFS (^,g,g,ff,g fe ) 
-2Re[Y qf j. -» Y-fj.~, L qk R qk ]m^m g V SFFFS {q 3 , q, q, g, q k ) 
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+(\ Y q N i q*Rq k \ 2 + \Ygfr t %Ls h \ 2 )mlVgwFs(<ij>Q> ( l>9,Qk) 

- 2 (\ Y qN t q* I 2 + \ Y qN iqj \ 2 ) Rc [ L q k R q k \ m q m 9 Y S FTF S% ^ ^ 9. 9, Qk) 

-2Re[F g ^.-.y^.-.^LiJm^.m^m § V sTWs (^, q, q, g, q k ) 

+( y q mg ] Y qm qk + Y ^ q Y qN^ m N, m <i Y Fsss^ }• (5.16) 

In eq. (5.15), / is summed over the symbols e, ji, t, u e , v lll v Tl u, d, c, s, t, b, with rif = 1 for leptons and 3 for 
quarks, and in eq. (5.16) the symbol q is summed over u, d, c, s, t, b. The indices j, k,n are each summed over 
the appropriate ranges (1,2 for squarks, sleptons, charginos, and charged Higgs scalars, 1,2,3,4 for neutralinos and 
neutral Higgs scalars, including the Goldstone bosons) wherever they appear. The masses and couplings appearing 

on the right-hand side are always running renormalized DR parameters. In all of the self-energy functions appearing 
in eqs. (5.15) and (5.16), the external momentum invariant is s — m\ with an infinitesimal positive imaginary part. 
Similarly, the pole masses for charginos in the MSSM can be written as 

M l ~ *cM dt = ml + JL-fig + p^ng, (5-17) 

for i = 1,2, where m 2 ^ are the tree-level squared mass eigenvalues, and the one-loop part is [61, 62], 

fig = n f (\Y fdiP , | 2 + \Y Wj \ 2 )B FS (f, Fj) + 2n f Re[Y fdtP ,Y Wj }m d m f B^ s (f, F 3 ) 

+ (\ Y CtNAl I' + \ Y C-N 3 4>t \ 2 )^F S {N 3 At) + 2 M Y ctN^- Y C-N j4 ,i] m C ^N^FS^j, <ft) 

+ d y c+cr o| 2 + \ Y c+c-*l\ 2 )?>Fs{C 3 ,4>l) + 2M Y ctc-4>l Y cjcr^m 6 m d B Ts {C ,^ k ) 

+e 2 m^[10 - 6WIJ + g 2 (\Off + |OjJ| 2 )B*v(iV,-, W) + 2g 2 Re[Of i Of*]m e m fr B FV (N j , W) 

+ (9 2 +9 ,2 )[(\0%\ 2 + \0^\ 2 )BMC J ,Z) + 2Re[0^0^]m d m d BMC^Z)], (5.18) 

and the two-loop part involving the strong interaction is 

ng - 8.9. 2 {^(lnc,Q« I 2 + l^c%Q 3 \ 2 )fi(Ci, q, Qj) + Tte\Ytf t Q.Ytf t Q j ]me t m q U(C i ,q, Qj) 
- 2Rc i Y q c l Q; Y ^c tqk R Q 3 L ^ M SFFSF(Q 3 , Q, q, q k ,g) 

+ RC [ Y qC t Q* Y QC t q' k L Qi L n + Y qC % Q 3 Y QC t q k R *Q , R q k ] m C t m 9 M S F F SF (Qj > Q > lAk , §) 

- 2Re [Y q Q, Q* Yqq, -» Lq, Rq k + Y-g, q. Yqq, - fc L*- k ] mg, m q M SFFSF (Qj,Q, q, q k ,g) 
-2Rc[Y qd .g.YZ d i9k LQ.Rq k ]m q m Q M s jr FsF (Q j , Q, q, q k ,~g) 

+2Re[Y q Q.Q*Y^£ i - k RQ.Rq k + Y*^QYQg.-»LQ.Lg k ]m q mgM SFFSF (Qj,Q, q, q k , g) 

+ Rc i Y qC t Q] Y QC,ql R Qi R Sk + Y qC t Q 3 Y QC t q k L *Q j L \Mc, m Q m Q m 9 M SFFSF(Qj ,Q,Q,Qk, 9) 

+(l y 9 c%Q* I 2 + \ Y qC,Q 3 \ 2 ) Y FSSF F (q, Q 3 ,Q 3 ,g, Q) + 2Re[F g(5 .Q»y- 6 .g.]m ( 5.m g V FsSFF (9, Q 3 ,Q 3 ,g, Q) 

- 2Re i( Y q C^ Y *c t Ql + Y ^QY q C t Q k ) L Q 3 R l k ] m Q m 9 Y FSS-FF^ Qj,Qk,9, Q) 

-2Re[Y qd ,Q,Y- e .Q k (Lq.Rq^ + RQ.L*Q k )]m 6 .m q m Q m- g V FSSFF (q, Q 3 ,Q k ,g, Q) 

+(\ Y q c i Q* L Qk\ 2 + \ Y q c i Q j R ^\ 2 ) V SFFFs{Q 3 ,q,q,g,qk) + 2Re[Y q ^Q,Y- dt Q 3 }m d m q V SFFFS (Q 3 ,q,q,g,q k ) 
-2Re[F g5 .Q,Fg 6 .Q.i 9 - fc J?|Jm e .m § V SFFFS (Q j ,g, q,g, q k ) 

+(\ Y q c i Q* R <ik I 2 + lYqGtQjLqk \ 2 ) m2 q Y sTffs(Qo . ^ Q,9,Qk) 
- 2 (\ Y q c % Q; I 2 + \ Y qC t Q 3 \ 2 ) R 4Lq k R*q k }m q m g V SFFTs (Q v q, q,g, q k ) 
-2Re[y^ i Q,r^Q j L^i?gJm 6t mgm 2 V slW5 (Qj, q, q,g, q k ) 
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+ (Y q c t Q ] Y^Q k + Y* di QY q * &iQl )m di m q Y Tsss (q, Qj,Q k , Q n )] }■ (5.19) 



In the one-loop part eq. (5.18), n/ = 1 for leptons and 
rif = 3 for quarks, and the symbols (f,F) are summed 
over the 12 ordered pairs: (e,v e ), {v e ,e), (fj,, v^), (f^,/i), 
(t,is t ), (v t ,t), (u,d), (d,u), (c,s), (s,c), (t,b), (b,t), 
while in the two-loop part eq. (5.19) the symbols (q, Q) 
are summed over the last 6 of these. The indices j,k,n 
are each summed over the appropriate ranges (1,2 for 
squarks, sleptons, charginos, and charged Higgs scalars, 
1, 2, 3, 4 for neutralinos and neutral Higgs scalars, includ- 
ing the Goldstone bosons) wherever they appear. In all 
of the self-energy functions appearing in eqs. (5.18) and 
(5.19), the external momentum invariant is s — m 2 ^ with 
an infinitesimal positive imaginary part. 

The numerical values of the two-loop neutralino and 
chargino pole masses are rather sensitive to the values of 
the model parameters. Most often, they can be expected 
to be at least several tenths of a per cent, but larger in 
some regions of parameter space. A study of the numer- 
ical significance of these results, and other contributions 
to the neutralino and chargino masses that are implicitly 
given in section IV, is deferred to future work. 

VI. OUTLOOK 

In this paper, I have presented results for radiative cor- 
rections to the self-energy functions and pole masses of 
fermions at two-loop order. The main specific motivation 
for this work is to allow future experimental data on su- 
perpartner masses to be connected to hypotheses for the 
mechanism of supersymmetry breaking. However, the 
strategy used is designed to be more flexible, with po- 
tential application to any semi-perturbative theory that 
may appear, anticipated or not, at the TeV scale. 

The application to the gluino mass may be particularly 
crucial, because of the relative strength of the SU(3)c 
gauge coupling, and the color octet representation of the 



gluino. When one extrapolates the soft-supersymmctry 
breaking parameters to very high energies [63]- [67] using 
two-loop or even three-loop [68] renormalization group 
equations, the gluino mass can also have a quite strong 
effect on the determination and running of other param- 
eters. It has been found that the combination of the 
Large Hadron Collider and a future Linear Collider may 
be able to pin down the gluino mass to 1% or so [1]. In 
that case, the two-loop corrections to the gluino mass 
will definitely be required. In general, the two-loop con- 
tributions for other fermions are parametrically smaller, 
but still worth including on the grounds that theoretical 
errors should be made negligible if at all practicable, in 
order to cleanly isolate the experimental implications of 
new data. 

For reasons of relative simplicity (and not principle), 
the calculations presented in this paper have neglected 
the masses of vector bosons in the two- loop part. Al- 
though this is a quite adequate approximation for many 
purposes, it can and should be improved on in future 
work. Also, the calculations have been presented here 
in a general form, and require specialization; this can be 
non-trivial for reasons related more to fatigue in writ- 
ing and typing than to conceptual difficulty. However, 
this specialization of general results seems well-suited to 
symbolic manipulation programs. In any case, the more 
lengthy results for specific contributions to the neutralino 
and chargino masses, for example, might be better placed 
in the innards of computer codes of the type described 
in [64]- [67] rather than explicitly on paper. 

I am grateful to Dave Robertson for valuable conver- 
sations and collaboration on the two-loop self-energy in- 
tegral computer program TSIL (ref. [43]) and Howard 
Haber and Herbi Dreiner for useful comments on section 
IIC. This work was supported by the National Science 
Foundation under Grant No. PHY-0140129. 
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